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1.  INTRODUCTION 

It  is  well  recognized  that  the  mean  sea  level,  as  a  surface  of  reference, 
varies  rapidly  'in  time  and  space1,  e.g.  see  Lennon  (1974).  The  variation 
of  mean  sea  level  from  the  geoid,  or  the  'quasi-stationary'  non-tidal 
sea-surface  topography,  e.g.  see  Mather  (1978),  is  not  very  well  known 
but  is  of  the  order  of  1  meter  r.m.s.  (root  mean  square).  The  vertical 
geodetic  networks  in  different  continents,  which  are  respectively  referenced 
to  mean  sea  level  determinations  at  coastal  tide  gages,  therefore  refer 
to  different  equiDotential  surfaces.  The  problems  of  unification  of  vertical 
datums  have  been  briefly  reviewed  by  Rapp  (1980). 

Colombo  (1980a)  proposed  that  a  'World  Vertical  Network'  may  be  defined 
by  'a  set  of  estimated  potential  differences  among  benchmarks  situated 
in  various  continents'  (ibid.,  p.2).  These  potential  differences  were 
postulated  to  be  independent  of  the  precise  definition  of  the  geoid,  or 
the  precise  knowledge  of  the  zero  harmonic  of  the  earth's  gravity  potential. 

The  principle  of  Colombo's  method  is  more  fully  described  in  Section  2, 
but  in  outline,  he  considered  the  potential  difference  between  two  bench¬ 
marks,  say  8MA  and  8MB,  each  in  a  different  continent,  by  computing  the 
normal  and  disturbing  potential  at  each  1 ocation  .subject  to  the  uncertainties 
in  the  knowledge  of  the  earth's  gravity  field  and  station  (benchmark)  locations. 
The  estimated  potential  difference,  AW(BMA.BMB),  could  be  made  more  precise 
through  redundant  determinations  by  computing  the  potential  differences 
between  some  other  pairs  of  benchmarks,  Pi  and  Qj  in  the  respective 
continents;  and  using  these  later  differences  to  estimate  AW(BMA.BMB) 
by  considering  potential  differences  AW(BMA,Pi)  and  AW(BMB ,Qj )  obtained 
through  gravity  and  leveling  observations. 

Colombo  considered  the  normal  gravitational  potential  described  by 
a  spherical  harmonic  expansion  to  degree  and  order  20,  with  standard  devia¬ 
tions  of  the  potential  coefficients  ('imperfect'  model)  as  determined 
by  Rapp  (1978),  and  also  with  zero  standard  deviation  ('perfect'  model). 

The  disturbing  potential  was  computed  from  gravity  anomalies  approximately 
uniformly  distributed  at  a  spacing  of  about  0.4°  around  each  benchmark 
ir  a  'cap'  of  spherical  radius  5°  (10°  was  also  considered  in  some  cases). 
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With  four  benchmarks  at  hypothetical  locations  in  Australia  and  five  bench¬ 
marks  in  North  America,  with  standard  deviation  of  relative  geocentric 
position  as  0.15  m  and  the  standard  deviation  of  potential  differences 
uW(BMA,P..),  iW(BMB,Qj)  as  0.1  \/T~ kgal  .m,  where  l  is  the  length  of  leveled 
connection  in  10 3  km,  Colombo  (1980a,  p.33)  estimated  the  •.to  do  n  in  ; :  . 
of  the  Australia  -  North  America  Vertical  datum  connection  as  about  0.3 
kgal.m  for  the  'imperfect'  model,  and  about  0.2  kgal.m  for  the  'perfect' 
model, i.e.  the  low  degree  potential  coefficients  to  degree  and  order  20 
beina  known  perfectly  without  any  error.  Colombo  (19P  extended  the 
results  to  consider  a  0.3  m  r.m.s.  unknown  shift  of  tl  arigin  of  the 
(benchmark)  coordinate  system  from  geocenter,  when  the  'Limated  standard 
deviation  of  the  vertical  datum  connection  rose  to  a  be  4  and  0.3  kgal.m 
respectively  for  the  imperfect  and  perfect  gravitation  ’  .odel .  A  slightly 
lower  value  was  obtained  for  hypothetically  located  benchmarks  for  the 
North  and  South  America  vertical  datum  connection. 

The  precise  coordinates  are  now  available  of  the  Lageos  Satellite 
Laser  Ranging  (SLR)  network  (Smith  et  al.,  1982)  in  support  of  the  NASA 
Geodynamics  program  (NASA,  1981,  p.6;  NASA,  1982,  p.7).  Improved  estimates 
are  also  now  available  of  the  potential  coefficients  in  the  spherical 
harmonic  expansion  of  the  earth’s  gravitational  field  (Rapp,  1981)  to 
degree  and  order  180.  The  present  study  applies  the  ideas  of  Colombo 
(1980a,  1980b)  to  the  currently  available  data,  which  is  described  in 
Section  3,  to  determine  the  current  accuracy  estimates  of  the  vertical 
datum  connections.  The  West  Europe-United  States  vertical  datum  connection 
is  examined  in  detail  in  Section  6,  as  a  greater  number  of  SLR  stations 
are  located  in  these  areas.  The  current  accuracy  estimates  are  also  deter¬ 
mined  for  some  other  vertical  datum  connections  in  Section  7. 

The  present  study  first  investigates  in  Section  4  the  optimum  require¬ 
ment  of  the  density  and  extent  of  gravity  anomalies  in  the  cap  around 
each  benchmark  for  estimating  the  disturbing  potential  at  the  cap  center. 
Anv  reduction  in  the  spherical  radius  of  the  cap  from  5'  makes  the  grav'T.v 
anomaly  data  acquisition  and  reduction  more  feasible.  This  reduction 
in  *;he  cap  radius  is  first  examined  in  terms  of  the  accuracy  of  ■astimat 


disturbing  potential.  Further  reduction  in  cap  radius  is  examined  in 
Section  6  in  terms  of  the  accuracy  of  the  vertical  datum  connection,  as 
a  larger  number  of  cap  pairs  may  be  available  for  a  reduced  cap  radius. 
Finally,  an  optimal  minimum  size  of  cap  radius  is  recommended  for  a  very 
modest  requirement  of  gravity  anomaly  data  while  still  ensuring  reason¬ 
able  accuracy  estimates. 

The  hypothetical  benchmark  locations  in  Colombo's  studies  (1980a, 
1980b)  were  placed  away  from  coastal  areas  so  that  gravity  anomalies  may 
not  be  required  over  marine  areas,  as  such  anomalies  would  be  of  lower 
precision  as  compared  to  those  obtainable  over  land  areas.  However,  as 
several  stations  of  the  SLR  network  are  located  in  coastal  areas,  the 
effect  of  lower  precision  of  gravity  anomalies  is  examined  in  Section 
4.5. 


The  accuracy  of  potential  difference  estimates  AW(BMA,Pi)  and  AW(BMB,Q 
obtained  through  leveling  could  be  of  some  concern  over  intra-continental 
distances,  e.g.,  see  Kumar  and  Soler  '1981);  Lachapelle  and  Gareau  (1980, 
Section  4).  The  effect  of  leveling  errors  on  the  vertical  datum  connection 
are  examined  in  Sections  6  and  7  of  this  report. 

The  effects  of  the  unknown  shift  of  the  origin  of  SLR  coordinate 
system  from  the  geocenter , and  the  random  position  errors  of  the  SLR  stations 
are  also  examined  in  Sections  6  and  7.  The  accuracy  of  the  normal  gravi¬ 
tational  field  is  considered  in  Section  4.4. 
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MATHEMATICAL  MODEL  FOR  THE  VERTICAL  DATUM  CONNECTION 


We  summarize  here  the  proposal  of  Colombo  (1980a,  1980b/  fur  es¬ 
tablishing  the  vertical  datum  connection  between  two  continents  or  regions 
A  and  B  by  estimating  the  potential  difference  AW(BMA,BMB)  between  two 
benchmarks  8M  A  and  BM  B  in  the  two  regions.  We  assume  that  there 
are  several  benchmarks  or  stations  and  Qj  respectively  in  the  two 
regions,  and  that  the  potential  differences  AW,(BMA,P.)  and  .\W,(BMB,Q.) 
have  been  obtained  through  leveling  (and  hence  the  subscript  i1  ).  Note 
that  P.  and  Q.  respectively  include  BMs.AandB,  when  obviously 
AW  ,  (BMA  ,BMA )  and  AW^BMB  ,BMB)  are  zero. 

2  . _1  Accuracy  Est i mate  of  the  Vert  i c_a  1  _  Da  t urn  Con n ec t i on 

Denoting  the  gravity  potential  W  as  a  sum  of  the  normal  gravi¬ 
tational  potential  U  ,  the  rotational  potential  •*  ,  and  the  anomalous 
potential  T  ;  and  also  introducing  errors  by  the  notation  ,  we  have 


AW(BMA,BMB)  -  U(Pi )  +  T(Pi)  +  i{P.)  +  AW£(BMA,P.) 

-  CU(Qj)  +  T(Qj)  +  0(QJ )  +  AWjfBMB.Qj)]  +  v  .. 


(2.1) 


where  v.j  =  *-U(Pi )  -  -UtQJ  +  tT^)  -  T(Qj)  +  -AWf  (BMA  ,P. ) ,  (2.2 

-  .AW,,  (BMB  ,Q j  ) , 

and  we  are  generally  following  Colombo's  (ibid.)  notation  for  ease  of 
reference.  The  errors  in  position  of  SLR  stations  propaaate  to  errors 
in  LI  and  I  ,  but  as  the  magnitude  of  errors  '  is  negligibly  small 
(•-  0.3.  )  as  compared  to  --U  ,  the  terms  s  have  been  omitted  in  (2.2' 


It  is  clear,by  a  simple  consideration,  that  the  maximum  number,  N  , 
of  linearly  independent  equations  (2.1)  between  paris  of  benchmarks 
P-  ,  Qj  for  establishing  the  vertical  datum  connection,  is  one  lev  than 
the  total  number  of  benchmarks  or  SLR  stations  considered  in  the  two  iciici 
"he  actual  choice  of  pairs  is  practically  determined,  thuuoh  a  -hitr.m  i !  . 
tor-  the  simplicity  of  design  (see  Section  6.1). 


Equation  (2.11  may  be  written  in  matrix  form  as: 


a  AW(BMA ,BMB)  =  £  +  v  (2.3) 

where  a,  £,  v  are  all  vectors  of  dimension  Ng.  a  is  the  design 
vector  with  all  components  equal  to  1,  AW(BMA.BMB)  is  the  unknown  potential 
difference  for  the  vertical  datum  connection  in  the  overdeterrnined  system 
(2.3);  £  is  the  vector  of  'observed'  potential  differences,  and  v  is 
the  vector  of  residuals. 

The  accuracy  estimate, c  AW(BMA,BMB),  of  the  vertical  datum  connection 
is  obtained  from: 

Ne  Ne 

c  AW(BMA,BMB)  =  (aT  V* 1  a)'1/2  =  [  Z  Z  (V-1),  A'l/2  (2.4) 

k  1  kl 

the  last  form  of  the  equation  being  a  consequence  of  all  components  of 
a  being  1.  In  (2.4),  V  is  the  variance-covariance  matrix  of  the  'ob¬ 
served'  potential  differences: 

V  =  V(cAU)  +  V(cAT)  +  V(sAAW  )  (2.5) 

where  the  right  hand  side  terms  are  the  variance-covariance  matrices  of 
the  errors  in  U(Pi)  -  U(Qj),  f(P.)  -  f(Qj),  AW£(BMA,Pi)  -  AW£(BMB,Qj) 
respectively.  The  computation  of  V(;AU)  and  V(eAAW£)  is  described  in 
Sections  2.2  and  2.3,  while  V(cAT)  is  considered  in  Section  5  after  con¬ 
sidering  the  accuracy  estimates  of  the  predicted  disturbing  potential 
in  Section  4.1.  The  notation  T(P^)  or  T(Q^)  indicates  that  'modified' 
gravity  anomalies  (Section  4.2)  have  been  used  for  predicting  T. 

It  is  also  clear  from  subsequent  discussion  that  the  errors  eAU,  £aT, 
eaaW^  are  primarily  due  to  different  parameters.  The  correlations  between 
these  errors  have  therefore  been  neglected  in  (2.5). 

2.2  Variance  Covariance  Matrix  of  Normal  Gravitational  Potential  Differences 

The  normal  gravi tational  potential  U  at  a  station  P  with  geo¬ 
centric  coordinates  (r,  s',  .* )  may  be  denoted  by: 


u(r, •  «tl.  I  (a/r-)n  7  P  (sin  i'JiC  cosna 
n=2  m=0 


+  S  sin  rn’k  '•  1 
nm 


(2.5) 


where  GM  is  the  gravitational  constant  times  the  mass  of  the  earth; 

N  is  the  highest  degree  of  the  spherical  harmonic  expansion;  a  here  is 
the  equatorial  radius,  ^nm»  ^nm»  ^nm  are  respectively  the  fully  normal¬ 
ized  Legendre's  functions  and  the  potential  coefficients  of  degree  n  and 
order  m. 

The  errors  -U  in  (2.6)  are  due  to  errors  in  station  coordinates, 

in  the  value  of  GM  ,  and  in  the  potential  coefficients,  i.e  .c  C  S  . 

nm  nm 

We  will  consider  the  effect  of  errors  ►: Cnm,  -Snm  entirely  while  evaluating 
the  error  T  in  Section  4.1;  and  consider  their  effect  on  cU  here 
as  zero.  The  effect  of  oGM/r  will  be  almost  zero  on  the  potential  dif¬ 
ferences  as  such  errors  cancel  in  the  differencing.  And,  as  already  mentioned, 
the  errors  of,  eX  for  the  SLR  stations  will  have  negligibly  small 
effect  on  eU  ,  as  compared  to  the  effect  of  or  on  o(J.  We  may  thus 
presently  consider  only  the  effect  of  errors  cr  ,  i.e., 

-r(P.)  _-.-ver(P.)  (2.7) 

where  it  is  sufficient  to  consider  an  average  value  v  of  normal  gravity 
on  the  ellipsoid.  This  approximation  is  reasonable  as  the  latitudinal 
spread  of  the  SLR  network  is  not  large  (see  Table  3.1).  For  simplicity, 
we  further  assume  for  this  study  that  the  standard  deviation  of  radial 
position  errors,  oer  ,  has  the  same  value  for  all  stations  of  the  SLR 
network.  For  actual  vertical  datum  connections,  it  is  not  necessary  to 
make  these  simplifying  assumptions. 

It  is  also  reasonable  to  assume  no  correlation  for  cr  in  the  deter¬ 
mination  of  the  positions  of  SLR  stations,  and  we  thus  have  a  diagonal 
matrix  for  V(fdsU)  with  each  diagonal  term  j2cAU  given  by: 

'•cLU  =  27-  r  .r  {[■..' 
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We  now  need  to  also  consider  the  additional  contribution  to  V(  AU) 
of  the  unknown  shift  of  the  SLR  coordinate  system  from  the  geocenter. 

The  error  er(P.)  caused  by  the  shift  (AX,  AY,  AZ)  of  the  origin  of 
the  coordinate  system  is  given  by  the  projection  of  the  shift  on  the  radius 
vector  r(P. ): 


X-A  X+  Y  •  AY  +  Zj AZ 

sr(Pi)  =  - - 1 - L- 

r(Pi) 


(2.9) 


where  (x^,  Y^ ,  zp  are  the  coordinates  of  P-. 


Denoting  the  coordinates  of  Q.  by  ( X  - ,  Y,,  Z.),  and  the  spherical 

s J  J  J  J 

distance  between  P.  and  Q.  by  '^(P.  ,Q^);  the  covariances  in  V(t;AU) 
in  addition  to  (2.8)  are  obtained  (see  Colombo,  1980b,  p.99,  equation 
(21))  by  using  (2.7)  and  (2.9): 

X.X.  +  Y.Y.  +  Z.Z. 

Y2  cov  ( £r(P  • ) ,  er(Qj)  =  Y2  - J A  '?2A 

1  J  r(p.)  r(Qj) 

=  Y2  ^  cos(^(Pi ,Qj))  (2.10) 


where  the  simplifying  assumptions  have  been  made  that  variances  of  shifts 
(g2AX,  a2AY,  a2 AZ)  along  coordinate  axes  are  not  correlated  with  each 
other,  and  are  all  equal  to  some  constant  a‘A  ; 

o2AX  =  o2 AY  =  o2 AZ  =  o2A  =  o2 s/3  (2.11) 


and  os  is  the  standard  deviation  of  the  shift  of  the  coordinate  system 
origin.  A  range  of  values  may  be  tried  for  os  ,  say  from  0  to  50  cm. 


2.3  Variance  Covariance  Matrix  of  Leveled  Potential  Differences 

The  standard  deviation,  c  ,  of  the  first  order,  class  I  aouble- 
run  leveling,  based  on  tolerances  (Federal  Geodetic  Control  Committee, 
1980,  pp.24  and  28)  between  forward  and  backward  measurements,  may  be 
quoted  as  0.5VTmm,  where  £  is  the  leveled  distance  hi  kilometer; 
(Vanicek  et  al . ,  1980,  p.  507).  However,  this  may  not  be  realistic  for 
long  lines  of  leveling  observed  over  a  period  of  several  months.  For 
the  determination  of  jW ^  over  intra-continental  distances  between  SLR 
stations,  a  more  conservative  estimate  (Personal  Communication  from 
C.T.  Whalen,  NGS  to  R.H.  Rapp,  June  1981)  would  be  o  =  2.  lO'VF  kgal  .m. 

Now,  if  H  is  expressed  in  103  km,  c=  .06VFkgal.m,  which  is  comparable 
to  0-lvT~kgal.m  used  by  Colombo  (1980a, b).  To  get  an  insight  into  the 
effect  of  leveling  errors  on  the  vertical  datum  connection,  a  range  of 
values  will  be  tried  for  o  =  cc£W,  in  Sections  6  and  7: 

X, 

=  k(0.ivf)  kgal  .m,  l  in  103  km,  k  =  0,1, 2, 4  (2.12) 

If  AM.  is  obtained  from  BMA  (or  BMB)  to  P.  (or  Q.)  by  separate 
K  „  1  J 

routes,  the  variances  of  would  be  uncorrelated.  The  leveling 

errors  covariance  matrix  in  one  region  would  have  off-diagonal  elements 

if  the  level  route  proceeds  to  a  station  through  a  preceding  station. 

For  example,  if  in  region  B,  there  are  five  stations  Q.  (including  BMB, 

which  is  numbered  as  1),  and  level  routes  are  independent  from  BMB  to 

each  station  except  to  station  5,  which  proceeds  through  station  4,  then 

V(-JiW?  )R  may  be  written  in  the  form  of  following  matrix: 


.01  k: 


V(-:AWf  ) 


B 


(1,1) 


(1,2) 


(1,3) 


(1.4) 

(1.4) 


/  1  A  ' 

V  i  ,4; 

( 1 , 4  ) +  ( 4 , 5 ) 


(2.13) 
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with  other  elements  zero.  The  term  (1,4),  for  example,  is  the  distance 
between  stations  1  and  4  in  103  km;  and  similarly,  other  non-zero  terms 
in  (2.13).  The  form  of  V(cAW^)  is  discussed  further  in  Section  6.1  (also 
see  Table  6.1). 

The  elements  of  V(eAAW^)  in  (2.5)  for  leveling  errors  in  the  potential 
difference  AW^(BMA.P^)  -  AW^(BMB,Qj)  would  be  formed  by  summation  of  cor¬ 
responding  elements  in  V(cAWp)A  and  V(eAW£)g  ,  according  to  the  selected 
pairing  of  stations  in  the  two  regions  in  the  N  equations  represented 
by  (2.1),  assuming  no  covariances  between  leveling  errors  in  region  A 
with  those  in  region  B. 
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3.  STATIONS  AND  GRAVITATIONAL  MODEL 


The  stations  used  for  computing  the  accuracy  estimates  for  the  vertical 
datum  connection  were  from  the  solution  SL4  (Smith  et  al.,  1982),  and  are 
described  in  Section  3.1.  The  gravitation  model  i:  described  in  ^  M'.i 
3.2.  Potential  differences  AW^  ,  through  leveling  and  gravity  observa¬ 
tions  in  a  region,  were  not  needed  (see  Section  2.3  for  leveling  error 
covariances)  for  the  vertical  datum  connection  accuracy  estimates.  But 
estimates  for  AW^  should  be  readily  available  when  actual  computations 
are  made.  A  very  wide  range  of  ocAW j  in  (2.12)  will  be  considered  for 
the  present  accuracy  studies  in  Sections  6  and  7. 

It  is  assumed  for  this  study  that  point  gravity  anomalies  will  be 
available  in  a  uniform  pattern  at  an  approximate  spacing  of  10'  (  18-20  km) 
around  the  SLR  stations  in  a  cap.  These  uniformly  patterned  anomalies 
may  be  predicted  from  existing  gravity  anomalies  in  the  area.  (See  Section 
4.1.  Additional  gravity  observations  may  be  required  in  some  areas). 
'Modified'  anomalies  (Colombo,  1980b,  Sec.  3.3)  are  envisaged  for  predicting 
T(P.)  or  T ( Q  • ) ,  which  will  be  discussed  in  Section  4.2.  The  choice  of 

*  J 

an  optimally  minimum  spherical  radius  of  the  cap  for  the  vertical  datum 
connection  would  be  discussed  in  Section  6.3. 

3.1  Lageos  Satellite  Laser  Ranging  Stations 

A  precise  set  of  coordinates  are  available  in  the  SL4  system  (Smith 
et  al.,  1982)  from  Lageos  laser  ranging  data  for  May  1976  till  August  1981. 

A  slightly  revised  list  of  coordinates  was  supplied  on  punched  cards  by 
D.  Christodoul idis  (Personal  communication,  June  1982).  Counting  only  as 
one  station,  when  different  nearby  sites  were  occupied  at  different  times, 
there  were  29  stations  in  the  solution.  Stations  7896  (Pasadena)  and  7110 
(Monument  Peak)  were  very  close,  about  1°  and  1/2 -  respectively,  to  stations 
7115  (Goldstone)  and  7062  (San  Diego).  After  omitting  7896  ano  7110.  the 
remaining  27  stations  have  been  plotted  in  Figure  3.1.  There  are  14  station 
jn  conterminous  United  States,  4  in  EuroDe,  2  each  in  Australia  and  South 
America,  2  each  in  the  Carribean  and  Pacific  Islands,  and  1  in  Hawaii. 
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Figure  3.1:  Lageos  Satellite  Laser  Ranging  Stations  in  the  SL4  System.  (Smith  et  al.,  1982). 
Stations  7110  and  7896  have  not  been  plotted. 
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_ker  Optical  Center;  X  =  Survey  Pad.  “  The  Geodetic  Coordinates  are  for  the  SL~  Solution 
,  1982),  with  Ellipsoid  Parameters;  a  -  6,378,144.11  in,  1/f  =  298.255 


Because  the  largest  number  of  SLR  stations  are  in  Europe  and  U.S., 
the  vertical  datum  connection  between  these  two  regions  has  been  examined 
in  detail  in  Section  6.  The  coordinates  of  these  18  stations  have  been 
listed  in  Table  3.1.  The  coordinates  refer  to  the  Laser  tracker  optical 
center  at  the  four  sites  in  Europe  and  the  SAO  site  at  Mt.  Hopkins,  Arizona. 
The  coordinates  at  the  remaining  11  sites  refer  to  the  survery  pad. 

3.2  Model  of  the  Earth's  Gravitational  Field 

Rapp  (1981)  computed  the  potential  coefficients  in  the  spherical  har¬ 
monic  expansion  (2.6)  of  the  earth's  gravitational  field  to  degree  and 
order  180.  The  expansion  was  in  fact  carried  out  to  degree  300,  but  there 
was  some  concern  (ibid.,  p.  24)  about  the  sharp  discontinuity  at  degree 
180  in  the  quadrature  weights  for  computing  the  potential  coefficients. 

The  data  sets  utililized  for  developing  this  field  were  an  a-priori  set 
of  potential  coefficients  to  degree  36  based  on  a  number  of  the  latest 
available  solutions  including  a  substantial  number  of  resonance  terms, 
a  l°xl°  anomaly  field  derived  from  the  Seasat  data  set,  and  the  latest 
available  l°x  1°  terrestrial  field.  In  principle,  the  three  data  sets  were 
used  to  provide  a  global  set  of  adjusted  l°x  1°  anomalies,  from  which  the 
final  set  of  potential  coefficients  were  obtained  by  the  quadrature  formulas. 

The  standard  deviation  (eCnm,  eSnm)  of  the  potential  coefficients 
was  estimated  considering  the  data  noise,  i.e.  the  r.m.s.  standard  deviation 
of  the  global  anomaly  set,  which  was  taken  as  20  mgals;  and  the  sampling 
error  due  to  the  finite  size  of  l°xl° blocks. 


The  effect  of  data  noise  is  several  times  larger  than  the  effect  of 
sampling  error  on  the  standard  deviation  of  the  potential  coefficients 
(Rapp,  1981,  p.  29).  Figure  11  (ibid.,  p.  33)  shows  the  coefficient  per¬ 
centage  error,  %En  ,  as  a  function  of  degree  n.  This  has  been  reproduced 
as  Figure  3.2  and  shows  that  %En  is  ahout  60%  at  degree  60  and  rises 
to  about  150%  at  degree  180. 


(eC  +e:S 
mn  nm 


)/(2n+l)]*/[  I  (rnni+^m)/(2n  +  l)]'-100 

m=0 


(3.1) 
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The  potential  coefficients  in  Rapp  (1981)  will  be  used  in  this  study 
to  describe  the  earth's  normal  gravitational  field.  The  coefficients  set 
to  degree  180,  or  a  lower  degree  subset,  will  be  termed  as  December  1981 
potential  coefficients  field  when  the  standard  deviations  as  given  in  ibid, 
are  used.  It  will  be  termed  as  ‘perfect  to  N  =  1C ,  20,  30'  if  the  standai d 
deviations  of  the  coefficients  to  degree  10,  20,  30  are  set  to  zero.  This 
corresponds  to  the  characterization  by  Colombo  (1980  a,b)  of  the  'imperfect' 
and  'perfect'  model  of  the  normal  gravitational  field. 


By  using  (2.6),  the  anomalous  potential  T(P)  at  a  station  P  ,  with 
geocentric  coordinates  (r,p',\),  may  be  expressed  by: 


T(P)  I(a/r)n  l.  Pr Jsin  7)  nn,c°sm''+>:$n  sinrm] 


N 


n 


r  n=2 


m=0 

n 


nm 


nm 


nm 


*  nXl  (a/r)nm»0  P"”(Sin  **  <3-2) 


where  we  have  now  considered  the  errors  >  -^nm  in  t^ie  norma^  gravi¬ 

tational  field  to  degree  N  in  defining  T(P)  . 

We  may  similarly  write  an  expression  for  the  gravity  anomaly  T.g(P) 
using  the  spherical  approximation  relation  (Heiskanen  and  Moritz,  p.  89): 


To  =  -  i  •  T 

nr  r 


(3.3) 


The  covariances  involving  T  and  Tg  at  points  P  and  P1  are 
then  given  by  (Colombo,  1980a,  p.  13): 


cov (T(P)  ,Tg(P ! ) )  = 


G  • 


N 


I  (  tt;  ) n  (2n+l)  (n-  1)  P  (cos  y  pp  , ) 


rr 


r.-2  rr 


n  n 


TP 


\i=fi+l(rr~)n  (2n^l)(n-l)  ?n(cos  vppl).  (3.4) 


(T(P)  ,T(P  1 ) )  =  (rrr)"(2n+l)  c;  P  (cos  -DD,) 


rr  „  rr 
n  =  2 


+  I  (—r)"  (2n+l)  ^  P n(cos  0 ppi)‘ 
n=N+l  rr  n  n  pk 


(3.5) 


cov(Ag(P),Ag(P'))  l  iyp-f  (2n+l)  (n-  1)  2  -:2  Pn(cos  yppi) 

n=2 


+  I  (£rr)n  (2n+l)(n-l)2  :2  P  (cos  tpp,)}  (3.6) 


n=N+l 


where 


;'  =  T  (c-2?  +  r2S  )/(2n  +  1) 

n  LnK-  nm  nm" v  ' 

m=0 


(3. ') 


and  the  unknown  c2C  .  e2S  may  be  approximated  by  the  estimated  variances 
:  2  nm  nm 

e(T  ,  eS  of  the  potential  coefficients,  c  is  the  r.m.s.  variation 
nm  nm  n 

of  the  potential  coefficient  by  degree: 


=„  -  C  l  (C=  *5’  )/(2->*  1)]* 


n  "mtn  nm  nm7 
m-0 


(3.8) 


which  is  related  to  anomaly  degree  variance  cn  by  (Jekeli,  1978,  p.  19): 


cn  =  y;  (n  -  l)2  (2n  +  1)  j2 


(3.9) 


The  anomaly  degree  variances  were  computed  by  the  'two-component1  global 
model,  using  the  spherical  approximation  *or  the  earth  with  radius 
R=r=r'  =6371  km: 


r n \  --  n-1  ,  :>n+2  ,  n- 1  ,  :Nn+2 

■(R)  "J1  ^  +  <2  Tn-2)  (nTpl  M 


(3.10 


with  the  Darameters  as  determined  by  Rapp  (1979,  d.  15,  Table  5,  Case  1' 


:  =  3.4050  mgal '  , 
fT 

i  =  —  *  0.998006, 
R3 


A  =  1  ,  B  = 


140.03  mgal " 

R3 

—  =  0.914232  (3.11) 

R: 

2 


The  advantage  of  modeling  the  degree  variances  by  equations  like  (3.10) 
is  that  it  allows  the  evaluation  of  infinite  series  in  the  second  part  of  the 
equations  (3.4)  to  (3.6)  by  using  finite  recursions  (see  Jekeli,  1978,  Sec. 
III. 4). 
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4.  ACCURACY  ESTIMATE  OF  PREDICTED  DISTURBING  POTENTIAL 


We  examine  in  this  section  the  accuracy,  -T  ,  of  predicted  T  at 
a  benchmark  (SLR  station)  using  the  gravity  anomalies  in  a  cap  of  radius 
v  .  The  variance-covariance  matrix  of  errors  in  T(P^)  -  T(Qj),  i.e. 

V(cAT)  in  (2.5),  will  be  discussed  in  Section  5.  The  evaluation  of  ;>;T 
is  described  in  Section  4.1,  and  the  'modified'  anomalies  ,*g*  used  for 
determining  ccT  are  described  in  Section  4.2.  The  effect  of  anomaly  density 
or  separation,  Av  ,  in  a  uniform  pattern  of  anomalies,  and  the  effect  of 
the  extent  or  spherical  distance  (radius)  f  of  the  anomalies  on  oLT  is 
examined  in  Section  4.3.  The  effect  of  the  errors  in  the  low  degree  coef¬ 
ficients,  and  the  effect  of  maximum  degree  of  the  normal  gravitational  field, 
on  ccT  is  examined  in  Section  4.4.  Finally  the  effect  of  the  errors, 
ccAg*  ,  in  modified  gravity  anomalies,  on  aeT  is  examined  in  Section  4.5. 

4.1  Disturbing  Potential  at  Cap  Centers 

We  may  predict  the  disturbing  potential  T(P. )  or  T(Q^)  at  a  cap 
center  or  station,  P..  or  ,  in  principle,  by  using  least  squares  col¬ 
location  (Moritz,  1980,  Sec.  14): 

T(Pi>  (4.1) 

where 

-  =  (=zz  +  =)_1-Tz  (4'2) 

and  the  estimated  accuracy,  :-T  ,  by: 

:2£T(P.)  *  Cn  -  fTCTZ  (4.3) 


r  -  r 

-TP i ,Ag  -  -Tz 
,g  with  T(P.) 

vector  consists 

dom  noise  part 


is  the  covariance  vector  of  the  gravity  anomalies  vector 
,  and  its  elements  are  given  by  (3.4).  The  Gravity  anomalies 
of  a  signal  part,  or  true  value  z  ,  and  uncorrelated  ran- 
n.  The  variance  covariance  matrix  of  the  signal  is  given 
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by  C22  =  (C  ^  );  and  that  of  the  noise  by  Q  ,  which  may  be  assumed  to 

be  a  diagonal  matrix,  and  computed  as  in  (4.5)  below.  Cyy  is  the  co- 
variance  (T(P  . )  ,T(P  . ) ) ,  which  was  assumed  to  be  the  same  at  all  benchmarks 
(cap  centers),  when  the  anomalies  have  the  same  uniform  pattern,  and  .  , 
for  each  cap.  CTT  and  elements  of  C  may  be  computed  by  (3.5)  and  (d. u) 
respectively.  The  computation  of  elements  of  Qzz  and  Cy2  are  modified 
as  in  (4.6),  (4.7)  below  when  using  ring  averages  of  gravity  anomalies. 

The  computations  of  covariances  Cy2  and  Qzz  are  done  much  faster 
if  the  anomalies  are  located  in  a  uniform  pattern  (Colombo,  1979),  and  the 
number  of  elements  are  much  reduced  if  ring  averages  Ag^  are  used  instead 
of  individual  anomalies  Ag^m  in  the  kth  ring  of  specified  radius 

around  the  cap  center: 


N.  ,  ~^km 

k  m=l 


d  1 

akk  =  rp-  f  2. 

N  ;  2  km 


and  for  simplicity,  we  assume  for  this  study  that  ckm  =  ;cig  for  all  k 
and  m.  This  simplifying  assumption  is  not  necessary  for  the  computations 
of  the  actual  vertical  datum  connections. 


Tne  covariances  involving  ring  averages  are  related  to  covariances 
involving  individual  anomalies  by: 

,  Nk  N1 


_ _  1  *  ' 

cov(Ag  ,Ag  )  =  -T--  VL  l  cov(Aq,  ,Ag 

K  1  Vl  m=l  n=l  km  ln 


:ov(T(P.),Agk)  cov(T(P.)f  AgJ 


rhere  were  an  integer  number  of  rings  say  N. 


such  that: 


in  a  cap  of  radius 


The  distribution  of  anomalies  in  the  rings  was  one  anomaly  at  the  cap  center 
(zeroth  ring);  6,  12  and  24  anomalies  in  the  first,  second,  and  third  rings. 
The  number  of  anomalies  thereafter  doubled  as  the  radius  of  the  ring  doubled 
(i.e.  in  the  6th,  12th,  18th,  etc.  ring),  and  staying  constant  otherwise 
(see  Colombo,  1980  a,  p.  23).  The  prediction  of  T(P.)  and  T(Qj)  was 
done  from  ring  averages  of  gravity  anomalies  by  utilizing  (4.5)  to  (4.7). 

It  is  assumed  for  this  study  that  the  'observed'  anomalies  at  separation 
of  nearly,  or  slightly  less  than,  ,  and  of  extent  slightly  larger  than 
(o  +  AcO  are  available,  though  these  are  not  located  in  the  uniform  pattern 
described  above;  e.g.  see  NOAA  (1982)  for  the  conterminous  United  States. 
Additional  gravity  observations  may  be  required  for  some  caps.  The  data 
requirement  is  very  modest,  as  will  be  discussed  in  Section  6.3.  This 
set  of  available  anomalies  may  then  be  used,  as  a  first  step,  to  predict 
anomalies  at  locations  giving  the  uniform  pattern  of  anomalies  described 
above.  We  now  discuss  in  Section  4.2  a  set  of  'modified'  anomalies  Ag* 
in  the  uniform  pattern;  the  modification,  in  principle,  lies  in  using  each 
cap  center  as  the  'leveling  datum'  for  all  anomalies  inside  that  cap.  This 
modification  bypasses  the  requirement  of  accurate  physical  realization 
of  the  geoid  using  coastal  tide  gages,  i.e.  the  problem  of  sea  surface 
topography. 

4.2  Modified  Gravity  Anomalies 

We  define  (Heiskanen  and  Moritz,  Sec.  8.3)  the  gravity  anomaly  -'.g 
at  a  point  P1  on  the  terrain  from  the  observed  gravity  g: 

-g(P')  *  g(P')  -  v(Q')  (4.8) 

where  7  is  the  normal  gravity  on  the  point  Q'  on  the  telluroid 
corresponding  to  P '  ,  such  that 

U(Q')  +  f(Q')  =  W(P')  (4.9) 

where  ‘(O')  t(P')  may  be  using  the  notation  of  (2.1).  However,  the 
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knowledge  of  W(P')  is  subject  to  errors  in  the  leveled  potential  dif¬ 


ferences  AW  (0,P')  =  W(0)  -  W(P')  over  long  distances  from  the  geoid 
potential  W(0)  determined  at  the  tide  gages,  with  additional  error  due 
to  the  unknown  sea  surface  topography. 

To  avoid  this  problem,  Colombo  (1980b,  Sec.  3.3)  proposed  the  use 
of  'modified1  gravity  anomaly  Ag*(P')  referenced  to  the  normal  gravity 
potential  through  the  corresponding  cap  center  P.  ,  such  that: 

Ag*(P')  =  g(P')  -  y(Q")  ( 


where  Q"  is  a  point  defined  on  a  geopotential  surface  with  the  normal 
gravity  potential 


U(Q")  +  4(Q")  =  U(Pi)  +  4-^)  -  AW  (P^P1) 


which  is  well  defined  if  we  observe  the  leveled  potential  differences 
Aw, (P.,P')  over  the  short  distances  (P^,P')  inside  a  cap  from  the  cor 
responding  cap  center.  We  thus  need  estimates  of  AW^(P..,P')  for  all 
gravity  anomalies  inside  a  cap  to  compute  the  modified  gravity  anomalies 
Ag*  in  (4.8)*.  Errors  in  AW^  would  propagate  into  errors  into  Ag*. 
These  would  of  course  be  much  smaller  than  the  corresponding  errors  in 
(4.8),  due  to  shorter  distances  and  will  also  be  free  of  the  errors  of 
sea  surface  topography.  The  last  statement  is  applicable  as  the  normal 
gravity  potential  at  the  cap  center,  Pi  ,  is  computed  from  the  SLR 
station  coordinates. 


We  further  note  that  as  the  right  hand  side  of  (4.9)*  equals 
W'P')  -  T(P.)  ,  and  assuming  T(Q")  =*  MQ'),  we  have  using  (4.9): 


U(Q')  =  U(Q" )  +  T ( P . ) 


Hence, 


A9*(P ' )  =  A9(P')  +  v(Q')  -  y(Q")  =  A9(P')  -  v* 

O'  f 

=  Ag(P')  +  T(Pi)  (4.11) 

as  a  spherical  approximation.  (4.11)  permits  a  correction  to  be  computed 
in  (4.12)  below  to  the  estimated  T(P.)  from  the  modified  anomalies,  because 
in  (4.11)  Ag*  may  be  considered  as  the  original  anomalies  with  a  bias 
term  kT(Pi-)  ,  k  =  2/r(P.).  If  Ag  is  in  mgals,  and  T  is  in  kgal.m, 
k  *  (2  x  10  ) / (6 . 371  x  106)  *0.3. 

In  (4.1)  to  (4.3),  the  covariances  are  based  on  Ag  through  (3.4) 
to  (3.6)  and  (3.10).  We  now  consider  the  effect  on  T  and  oeT  >  if  the 
data  consists  of  the  modified  anomalies  Ag*  instead  of  Ag  ,  but  the 
covariances,  and  hence  f  in  (4.2),  is  based  on  Ag.  Let  the  estimate 
of  T(P.)  be  denoted  by  T(P.. )  or  T(P.)  based  on  the  data  Ag  or 
Ag*  respectively,  where  the  overbar  indicates  that  ring  averages 
are  used  as  in  (4.4)  and  (4.7)*,  Nq  =  Np  +  1.  Then: 

Ng  NP  Nq 

T(P  )  =  T(P  )  +  eT(P  )  =  l  f I Ag,  =  l  f, Ag*,  -  kT(P  )  l  f, , 

1  1  1  1=1  1  1  1=1  1  1  1  1=1 
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Henceforward,  whenever  we  use  T(P^)  or  jcT(P^)  ,  it  will  be 
understood  to  mean  T(P^)  or  c£T  ( P  ^ )  based  on  modified  gravity  anomalies 
Ag*  ,  but  all  covariance  would  be  based  on  Ag  through  the  anomaly  degree 
variance  model  (3.10)  used  for  evaluating  (3.4)  to  (3.6).  The  factor 
(1  +  kzf i )  and  (1  +  kzf^)2  is  then  required  to  be  used  as  in  (4.13) 
and  (4.13)  respectively. 


The  implementation  of  the  above  equations  was  done  by  slight  modifi¬ 
cations  of  the  programs  developed  by  Colombo  (1980a,  Appendix  B). 

4.3  Extent  and  Density  of  Gravity  Anomalies 

We  now  examine  the  effect  of  i.e.  the  cap  radius  or  extent  of 
gravity  anomalies  around  the  cap  center;  and  the  effect  of  Av  ,  i.e.  the 
density  or  separation  of  anomalies,  on  the  accuracy  estimate  c£T.  Colombo 
(1980  a,b)  had  considered  ^  =  5°  ,  but  any  reduction  in  i p  would  make 
the  realization  of  Ag*(P‘)  ,  with  the  needed  estimates  of  AW^(Pi,P') 
in  (4.9)*,  easier  in  practice. 

Figure  4.1  shows  a  plot  of  the  covariance  Cy  ^  =  cov(T(P),Ag(P')) 
in  (3.4)  for  variation  in  the  spherical  radius  ip  up  to  5°,  where  the 
normal  gravitational  field  is  represented  by  the  Dec.  '81  P.C.  field  ('im¬ 
perfect'  model)  and  also  the  'perfect'  model  to  N  =  30  (see  Section  3.2). 
The  curves  for  coefficients  perfect  to  N  =  10,20  were  not  plotted  to  avoid 
congestion;  they  fall  between  the  two  plotted  curves.  The  value  of  N 
in  the  expansion  of  T  in  (3.2)  and  the  corresponding  expansion  of  Ag 
in  (3.3)  was  taken  as  60.  The  use  of  this  specific  value  of  N  will  be 
discussed  in  Section  4.4.  The  plot  of  the  covariance  function  CT  .  woul 
of  course  be  more  damped  out  if  a  still  higher  degree  espansion  is  used 
in  (3.2)  and  (3.3),  i.e.  when  N  is  taken  as  90,  120,  150,  or  180.  It 
is  clear  from  Figure  4.1  that  CT  ^  decreases  rapidly  with  increase  in 
.  ,  and  the  extent  of  anomalies,  for  predicting  T  ,  may  be  reduced  con¬ 
siderably  from  j  =  5’  used  by  Colombo  (1980  a,b). 

For  a  specific  value  of  O  ,  we  have  also  to  consider  the  effect  on 
•  T  for  different  density  or  separation,  A A  ,  of  the  anomalies  (See 


COVARIANCE  CT>Ag=  [T  (P),  Ag  (P')]  IN  MGALS.2  M.IO 


SI® 


Figure  4.1:  Variation  of  Covariance  ( T , Ag )  with  Spherical 
Distance  <l>. 

December  1981  Potential  Coefficients  Field 
with  specified  standard  deviations. 

Perfect  to  Degree  30  means  standard  deviation 
of  Potential  Coefficients  up  to  degree  and 
order  30  set  equal  to  zero. 


DEC.  '81  POTL.  COEF.  FIELD 


PERFECT  TO  DEGREE  30 


5l  I  2  3  4  5 

SPHERICAL  DISTANCE  ^(P,P')IN  DEGREES 
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equations  (4.13)  and  (4.7*)).  The  values  of  '-:T  are  given  in  Table  4.1 
for  Ac  =  1/3°,  1/4°,  1/6°  for  some  values  of  >  .  The  total  number  of 
anomalies  needed  for  some  specific  combinations  of  Ay  and  y  are  also 
given.  The  normal  gravitational  field  is  described  by  Dec.  1 8 1  P.C.  field 
to  degree  N  <  60;  and  osAg*  is  2  mgals.  (The  effect  of  variation  in 
N  and  oeAg*  would  be  considered  in  Sections  4.4  and  4.5  respectively) . 


Table  4.1  Accuracy  Estimate,  ;?T  (Kqal.m),  of  Disturbing  Potential 
at  Cap  Center  for  Various  Density,  Ay°  ,  and  Extent,  .  , 
of  Uniform  Pattern  of  Modified  Gravity  Anomalies,  Ag*. 
jcAg*  =  2mgals;  December  '81  Potential  Coefficients  Field. 
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Table  4.1  confirms  the  obvious  expectation  of  lower  ocT  for  a  specific 
density  Ay  ,  as  a  is  increased;  and  also  of  lower  ccT  for  a  specific 
,  as  the  density  Ay  is  increased.  However,  a  large  y  along  with 
a  large  Ay  is  not  suitable  for  realization  in  practice,  as  the  total 
number  of  anomalies  in  a  cap  becomes  very  large,  e.g.  Ay  =  1/6°,  y  =4  or  5 

We  can  also  see  from  Table  4.1  the  possibility  of  an  optimum  com¬ 
bination  of  v  and  Ay  .  For  example,  v  =  3°,  Ay  =  1/6 "  gives  about 
the  same  value  of  cC T  .as  y  =  4°,  Ay  =  1/4°.  Again,  y  =  2°,  A  =  1/6 ' 
gives  slightly  higher  ceT  than  y  =  3°,  Ay  =  1/4 but  in  fact  lower  than 
A.  =  1/3°  for  either  y  =  4°  or  5C. 


It  is  obviously  preferable  to  have  a  smaller  v  ,  as  we  need  the  es¬ 
timates  of  AW  (P.,  P')  for  all  anomalies  in  a  data  cap,  and  these  can 


be  obtained  with  greater  precision  for  smaller  ..  .  Additional  tests  were 
carried  out  to  determine  ceT  for  ;  <  3°  ,  but  with  greater  density  of 
anomalies,  i.e.  Ay  _<  1/6°.  The  results  are  given  in  Table  4.2. 

Table  4.2  Accuracy  Estimate,  cc'T  (kgal.m), 

of  Disturbing  Potential  at  Cap  Center. 

Extent  ip  <  3°;  Ay  <  1/6°; 

Other  particulars  as  in  Table  4.1. 
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It  is  thus  adequate  to  have  the  anomalies  in  a  uniform  pattern  (4.7)* 
at  a  separation  Ay  ~  10 1 >  as  a  more  densely  spaced  set  does  not  lead  to 
significantly  lower  aeT  .  We  will  examine  the  accuracy  of  vertical  datum 
connection,  specifically  of  Europe  -  U.S.  connection,  using  y  of  3°,  2" 
and  lc  with  a  larger  number  of  equations  (2.3)  being  available  as  c  Is 
decreased. 

The  results  in  Table  4.1  and  4.2,  and  a  few  additional  results,  have 
been  plotted  in  Figure  4.2. 

4.4  Accuracy  and  Highest  Degree  of  Normal  Gravitational  Fi el d 

It  is  of  interest  to  know  how  much  the  accuracy  estimate  c-;T  would 
be  improved  if  the  low  degree  coefficients  in  the  gravitational  field 
were  known  to  a  higher  accuracy  than  at  present  (Section  3.2).  •• T  is 

accordingly  tabulated  in  Table  4.3  for  the  current  December  '81  field, 
as  well  as  for  the  cases  if  the  standard  deviation  of  the  coefficients 
is  zero  (perfect)  to  degree  N  =  10,  20  or  30.  The  extent  of  anomalies 
in  a  data  cap,  i.e.  :  ,  is  varied;  but  as  determined  in  Section  4.3, 
the  anomaly  separation  Ac  in  all  cases  is  taken  as  Av  =  1/6  . 


ACCURACY  o-rr  kgal. 


Figure  4.2: 


Variation  of  Accuracy  Estimate  o£T  of  Disturbing  Potential 
with  Anomaly  Cap  Radius  w  and  Anomaly  Density  aw. 

oeag*  =  2  mgals.  December  1981  Potential  Coefficients 
Field. 

The  Total  number  of  anomalies  in  a  data  cap  are  given 
in  Tables  4. 1  and  4.2. 
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Table  4.3  Accuracy  Estimate,  ::T  (kgal.m) 

of  Disturbing  Potential  at  Cap  Center. 
Variation  Due  to  Accuracy  of  Low 
Degree  Normal  Gravitational  Field. 

Anomaly  Separation,  Lv  =  1/6  .  ’oig*  =  2mgals 


Accuracy  of  Low  Degree  Normal  Field  | _  r-:T  (kgal.m) 
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The  results  in  Table  4.3,  and  some  additional  results  for  0  >  3C,  are 
plotted  in  Figure  4.3.  We  clearly  see  that  for  a  given  accuracy  of  the 
normal  gravitational  field,  jeT  improves  as  ^  is  increased.  However, 
oeT  improves  more  strikingly  for  a  given  0  ,  as  the  accuracy  of  low 
degree  coefficients  improves.  We  do  not  expect  a  perfect  gravitational 
model  in  the  near  future,  and  all  vertical  datum  connections  in  this  report 
in  Sections  6  and  7  are  computed  with  the  present  accuracy  of  coefficients 
in  the  December  '81  field.  However,  we  may  draw  two  conclusions  from 
Figure  4.3.  Firstly,  the  extent  v  of  ^g*  may  be  reduced  as  the  accuracy 
of  low  degree  normal  field  improves;  e.g.,  ~eT  of  about  0.4  kgal.m  may 
be  achieved  with  anomalies  in  a  data  cap  of  33,  2 3  or  1 with  Dec.  '81 
field,  and  a  field  perfect  to  N  =  10,  20  respectively .  Secondly,  we  may 
expect  limiting  values  of  cCT  of  about  0.3  kgal.m  with  data  cap  of 
-  =  1°,  which  will  improve  to  about  0.2  kgal.m  with  .  =  2C.  These  improvements 
will  come  primarily  from  better  knowledge  of  low  degree  coefficients, 
larger  data  caps  (v  >  3J)  or  greater  data  density  <  1/6  )  will  not 
contribute  sianif icantly  to  lowering  the  accuracy  ('-T)  of  disturbinn 
ootential  at  cao  centers. 
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Figure  4.3:  Variation  of  Accuracy  Estimate  ceT  of  Disturbing  Potential 
with  Accuracy  of  Low  Degree  Normal  Gravitational  Field. 
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It  is  also  of  interest  to  examine  the  effect  on  ■  T  in  (4.13)  of 
the  maximum  degree  N  of  normal  gravitational  field  used  in  defininc 
T  in  (3.2  );  and  the  covariances  (3.4)  to  (3.6)  required  in  comoutinq 
.-'T  in  (4.3).  This  is  shown  in  Figure  4.4  for  the  December  '81  coefficients 
field,  for  anomaly  separation  at  av  =  1/6 J,  v.tg*  =  2  maals;  and  for  data 
cap  radius  j  =  0.5°,  1%  2°  and  3°.  There  is  no  significant  difference 
in  :/T  in  any  of  the  cases  considered,  when  the  maximum  degree  N  of 
normal  gravitational  field  is  increased  beyond  60.  This  is  due  to  the 
fact  that  ther.m.s.  accuracy  en  and  the  r.m.s.  variation  of  the  coefficient 
•  (defined  in  (3.7)  and  (3.8)),  tend  to  be  of  the  same  order  of  magnitude 
when  N  60  (see  Table  3.2).  However,  there  is  a  marked  difference  mi 
the  computing  time  in  the  covariances  (3.4)  to  (3.6)  as  the  maximum  decree 
N  is  increased.  The  second  term  in  equations  (3.4)  to  (3.6)  is  ubtained 
by  finite  recursions  with  the  use  of  the  degree  variance  model  (see  end 
of  Section  3.2);  but  the  first  term  in  equations  (3.4)  to  (3.6)  requires 
the  summation  of  a  series,  which  becomes  expensive  in  computer  time  as 
N  is  increased.  For  ..  =  3°,  with  18  anomaly  rinos  (4.4),  the  computing 
time  for  :-:T  for  one  data  cap  on  the  Amdahl  V/6  computer  was  13.4  seconds 
for  N  =  60,  and  36.1  seconds  for  N  =  180.  However,  as  the  number  of  anomaly 
rings  are  reduced  to  12  for  i  =  2°,  the  computing  time  for  jlT  was  8.7 
seconds  for  N  =  60,  and  11.4  seconds  for  N  =  180.  The  difference  in  computinn 
time  is  further  reduced  when  <  2  . 

4.5  Accuracy  cf  Gravity  Anomalies 

The  accuracy  estimate  of  modified  gravity  anomalies  propagates  into 

■  T  through  (4.5)  which  defines  the  noise  matrix  used  in  computing  f 

in  (4.2)  and  •--cT  in  (4.3);  and  finally  through  (4.13)  we  also  take  into 

account  the  underlying  covariance  model  used  for  predicting  the  disturbing 

potential  T(P.)  or  i  ( Q  j )  at  cap  centers  in  (4.12).  This  propagation 

to  'IT  strictly  refers  to  the  random  errors  in  gravity  anomalies,  which 

deoends  on  the  precision  of  gravity  observations  and  the  position,  primarily 

the  height,  of  the  observation  point.  The  computation  of  modified  anomalies 

requires  TW.(P.,P‘)  observations  through  leveling  in  (4.9)*  and  (4.3)* 

from  the  cap  center  over  distances  of  the  order  of  .  =  3  ,  or  smaller. 

Hence,  a  conservative  estimate  of  2  meals  was  "sed  for  ■  'a*  *  for 

v  km 
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Figure  4.4:  Variation  of  Accuracy  Estimate  of! 

of  Disturbing  Potential  with  Highest 
Degree  of  Normal  Gravitational  Field. 


oeAg*  =  2  mgals.  Anomaly  Density 
All;  =  1/6°. 

December  1981  Potential  Coefficients 
Field. 


all  k  and  m  in  (4.5),  i.e.  for  simplicity  in  this  study,  all  gravity 
anomalies  were  assumed  to  have  a  standard  deviation  of  2  mgals.  This 
estimate  takes  into  account  the  gravity  observation  errors,  the  additional 
errors  in  leveling  from  the  cap  center  and  in  the  prediction  of  the  uniform 
pattern  of  anomalies  described  below  (4.7)*. 


However,  to  examine  the  effect  on  --T  of  larger  ccAg*  ,  several 
values  of  oeAg*  were  tried  up  to  10  mgals  for  different  cap  sizes.  These 
are  tabulated  in  Table  4.4. 


Table  4.4  Accuracy  Estimate,  'CT  (kgal.m), 

of  Disturbing  Potential  at  Cap  Center. 

Variation  due  to  Accuracy  Estimate,  ocAg*, 
of  Modified  Gravity  Anomalies. 

Uniform  Pattern  of  Anomalies  with  Density  Ac  =  1/6'. 
December  ’81  Potential  Coefficient  Field. 
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The  results  in  Table  4.4,  and  a  few  additional  results  for  the  case 
if  the  low  degree  (N  <  10)  potential  coefficients  were  known  perfectly, 
have  been  Dotted  in  Figure  4.5.  It  is  clear  from  Table  4.4  and  Figure 
4.5  that  *_T  depends  strongly  on  the  extent  of  gravity  anomalies  (with 
optimm  density  a.  =  1/6  );  and  T  is  not  particularly  affected  by  the 
variation  in  •  q*.  This  result  is  very  helpful  in  answering  the  concern 

about  the  aci.ura.-y  rtiaui remen t  of  gravity  anomalies  over  marine  areas 
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Figure  4.5:  Variation  of  Accuracy  Estimate  ofT  of  Disturbing  Potential 
with  Accuracy  of  Gravity  Anomalies. 
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around  those  SLR  stations,  which  are  located  in  coastal  areas  (see  F inure 
3.1).  Even  if  the  random  errors  of  some  modified  gravity  anomalies  are 
quite  large,  the  accuracy  estimate  ■■ T  of  disturbing  potential  at  cap 
center  is  obtained  to  about  '41,  '51,  '67  kaal.m  for  .  =3,2',  1  re¬ 
spectively,  if  modified  gravity  anomalies  are  available  in  a  uniform  pat¬ 
tern  at  density  Ay  =  1/6°. 

A  systematic  error  in  gravity  anomalies  due  to  height  error  of  as 
large  as  50  cm  is  shown  by  Colombo  (1980,  p.  20,  footnote)  to  cause  an 
error  of  the  order  of  0'05  kgal .m  in  2T  ,  which  may  be  neglected.  A 
constant  height  error  in  all  stations  will  result  in  a  constant  bias  in 
all  gravity  anomalies.  As  the  pattern  of  gravity  anomalies,  i.e.  the 
extent  v°  and  density  ,  is  the  same  around  each  cap  center,  the 

resulting  cT  will  be  nearly  the  same  and  will  largely  cancel  out  in 
disturbing  potential  differences  between  cap  pairs.  For  similar  reasons, 
a  systematic  error  in  GM  will  largely  cancel  out  from  the  normal  and  dis¬ 
turbing  potential  differences  between  cap  pairs.  The  effect  of  random 
errors  in  positions  of  cap  centers,  and  the  geocentricity  of  the  coordinate 
system  ,  on  :-:.£U  was  discussed  in  Section  2.2. 

We  now  first  examine  in  Section  5  the  variance-covariance  matrix 
of  disturbing  potential  differences  V('dT)  in  (2.5)  between  cap  pairs, 

P.j,  Q.  in  regions  A  and  B  .  The  optimum  extent  y  of  gravity  anomalies 
around  each  cap  center  is  then  examined  in  Section  6  by  T,'.W(BMA,BMB)  in 
(2.4)  for  the  case  of  Europe  -  U.S.A.  vertical  datum  connection,  in  view 
of  larger  number  of  equations  being  available  as  v  is  decreased  below 


5. 


VARIANCE  COVARIANCE  MATRIX  OF  DISTURBING  POTENTIAL  DIFFERENCES 


The  variance  covariance  matrix  V(Id.T)  in  (2.5)  for  errors  in 
i(P.)  -  l(Qj)  are  given  by  propagation  of  covariances  through  a ' -1 ( P . ) 
and  ' '  :2T(Q  - )  in  (4.13)  and  by  cov(eT(P.),  LT(Q-)).  For  example,  the 

J  1  J 

diagonal  element  vkk  of  V(cAT)  is  given  by: 


»kk  *  "'<?T(Pik>  -  n(Pjk))0 

-  ^T<Pik>  +  -  2M{rr<Plk>  'T(QJk)i. 


(5.1) 


where  M  is  the  averaging  operator,  and  k  =  1,...,  Ng  (see  equation  2.3)). 
Similarly,  the  off  diagonal  element  v ^  of  V(e£T)  is  given  by 
vk]  =  M((2T(Pik)  -  ?T(Qjk))  (?T(Pi])  -  CTCQ^))} 

=  Mi--T(P.k)  tJ{Pu)}  +  M(2T(Qjk)  CT(qjVi  -  M{CT(P-k)  CT(Q^); 

-  MuT(P11)  CT(Qjk)}(5.2) 


and  a  representative  term  M{-5T(P-k)  -T(Q^);  -  Mi  T(P.)  -:T(Qj.}  for 
convenience  of  notation,  and  using  (4.12),  is  given  oy: 


M  T(P.)  T(Q, 

*  J 


fT(fV  -  fI 

-  y  ■' 


[!t  F(ftl  ;f> 


.0.  T (0  . )  -  f'.a. 
•  j  -  r  j  j 


r(QjT  J 


-T 


^CTT(^!.0j)-2C-TJ(?i.0j)f  *  (H 


h  fi )"'1  (1+  rib  fjr' 

'  J 


where  J  are  the  vectors  of  rmo  averages  ut  modified  gravity  anomalies: 
'ner  notation-,  are  similar  to  (4.1)  to  (4.3);  and  we  have  assumed  in 
(0.3 '  fhar  ‘h*»  iat*-m  ,  f  a  noma  1  ,  ana  -.heir  noise  matrix  ,  is  the  same 

i  1,0  1  _..v  .);>  >.-n  f  <?'- . 


The  computation  of  covariances  QcLcL  of  ring  averages  of  aravity 
anomalies  in  (5.3)  is  very  expensive  on  computer  time,  if  these  are  com¬ 
puted  rigorously  by  (4.6)  and  (3-6).  It  is  adequate  for  (5.3)  to  compute 


the  covariance  of  ring  averages  k  and  1  ,  i  .e.  M  'ig^AQ^ } 
(1980,  p.  31,  (4.17)): 


by  Colombo 


M{Agik  Ag^}  =  M{^7 


HBiT  !  m(-V  '■>d'  •SaTli-  '  *<v,.  ••)<* 

'k  o  1  o 


T-  l-l  (2n+l)(n-l)2  -  P„(cos  %)  P„(-os  ,,)  P„(cos  flj) 


+  ni+1  (2n+1)(n'1):  "'n  Pn(cos  *k>  Pn<cos  v, )  P„(cos  (5.4) 


where  are  the  radii  for  the  rings  k 
spherical  distance  between  cap  centers  i  and 
are  as  in  (3.6) . 


and  1  ,  v.j  is  the 
j  ,  and  other  notations 


The  covariances  for  disturbing  potential  at  cap  centers  from  (5.3) 

have  been  plotted  in  Figure  5.1  against  spherical  distance  ecjj  when 

all  caps  are  of  a  given  size,  i.e.  3°,  2°,  1°,  0.5°.  December  1981 

Potential  coefficients  field  was  used  with  N  =  60  and  N  =  500.  The 

max 

elements  of  the  matrix  V(£AT)  for  disturbing  potential  differences 
are  computed  according  to  (5.1)  and  (5.2)  based  on  the  cap  pairs  P. , 
for  the  vertical  datum  connection  equations  (2.1). 


6.  ACCURACY  ESTIMATE  OF  EUROPE- USA  VERTICAL  DATUM  CONNECTION 


We  note  from  Figure  3.1  that  the  maximum  number  of  satellite  laser 
ranging  (SLR)  stations,  whose  coordinates  are  available  in  the  SL4  co¬ 
ordinate  system  (Table  3.1)  and  which  are,  or  may  be,  interconnected 
by  leveling  in  a  region,  are  14  in  USA  and  4  in  Europe.  We  therefore 
examine  in  this  section  the  various  possibilities  of  achieving  Europe- 
USA  vertical  datum  connection  by  estimating  the  disturbing  potential 
from  modified  anomalies  in  various  capsizes  around  the  SLR  stations.  The 
effect  of  errors  in  the  geocentric i ty  of  the  SL4  coordinate  system,  the 
effect  of  random  errors  in  SLR  station  positions,  and  the  effect  of 
leveling  errors  inside  each  region  on  the  accuracy  estimate  of  the  ver¬ 
tical  datum  connection  is  examined  in  Section  6.2.  We  then  discuss  in 
Section  6.3  the  optimum  capsize  for  Europe-USA  vertical  datum  connection. 

6.1  SLR  Stations  for  Vertical  Datum  Connection 

Figures  6.1  and  6.2  show  the  SLR  stations  available  in  Europe 

and  USA  for  the  vertical  datum  connection.  A  cap  of  radius  T  =  1°  has 
been  drawn  around  each  cap  center,  in  which  '  modified'  gravity  anomalies 
Ag*  (Section  4.2)  may  be  used,  at  a  spacing  of  10'  (AC  =  1/6°)  in  a  uniform 
pattern  (see  below  (4.7)*),  to  predict  T C P ^ )  or  T ( Q ^ ) .  With  c  =  1°, 
we  have  4  'benchmarks'  P^  in  Europe  and  14  'benchmarks'  Qj  in  USA  to 
establish  the  vertical  datum  connection,  by  using  Ng  =  17  independent 
equations  (2.1).  To  ensure  that  these  equations  are  linearly  independent 
(see  (2.3)),  and  to  ensure  numerical  stability  in  the  solution  of 
AW(BMA,BMB),  the  caps  should  not  be  overlapping  so  that  different  sets  of 
gravity  anomalies  are  used  to  predict  T  at  each  cap  center. 

We  recall  from  line  1  of  Table  4.2  that  the  accuracy  estimate  -*fT 
of  the  disturbing  potential  at  cap  center  improves  as  the  cap  radius 
is  increased.  However,  the  number  of  caps  then  becomes  smaller  to  avoid 
the  overlapping  of  caps.  Figures  6.3  and  6.4  show  the  available  cap 
centers  in  USA  when  the  radius  v  is  2°  and  3°  respectively.  The  avail¬ 
able  cap  centers  in  USA  are  therefore  reduced  from  14  to  10  and  5,  as 
is  increased  from  1  to  2°  and  3  .  All  the  4  caps  in  Europe  remain 
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Satellite  Laser  Ranging  Stations  in  USA  in  the  SL4  System  (Smith  et  a  1 . ,  1982). 

Gravity  Anomaly  Ag*  Cap  Radius  ip  =  1°  for  Predicting  Disturbing  Potential  at  Cap  Centers 


Cap  Centers  in  USA  for  Vertical  Datum  Connection.  Cap  Radius 


Cap  Centers  in  USA  for  Vertical  Datum  Connection.  Cap  Radius 


non-overlapping  even  with  y  =  3°,  except  for  a  very  slight  overlap 
between  the  caps  around  WETZEL  and  KOOLAS.  The  number  of  equations 
(2.1)  for  the  vertical  datum  connection  are  thus  =  17,  13,  8  in 
tne  three  cases.  The  advantage  of  increased  accuracy  in  the  predicted 
T(P^)  or  T ( Qj )  is  therefore  somewhat  offset  by  the  reduction  in  tne 
number  of  equations  N  in  (2.3).  Also,  the  increase  in  ,  very 
sharply  increases  the  number  of  modified  anomalies  (Table  4.2,  line  1, 
last  4  columns)  required  in  each  cap;  and  accordingly  increases  the 
field  work  for  determining  the  potential  differences  *iW;  ( P  ,  P ' )  in 
(4.9)*  from  the  cap  centers  to  the  anomalies. 

The  number  of  linearly  independent  equations  (2.1)  are  one  less  than 
the  tGtal  number  of  caps  in  both  regions.  The  choice  of  cap  pairs  P,, 

Qj  in  (2.1)  is  arbitrary  so  long  as  these  equations  are  independent; 
the  actual  choice  is  best  made  for  the  simplicity  of  vertical  datum 
connection.  For  example,  when  3  =  3°,  with  4  caps  P^  in  Europe  and  5 
caps  Qj  in  USA,  the  cap  pairs  for  forming  equations  (2.1)  for  the  vertica 
datum  connection  may  be  as  follows  (see  Table  3.1  for  details  of  stations) 

Pr-  WETZEL  GRASSE  KOOLAS  FINLAS 


Oj  ••  PLTTVX  OUINCX  FTDAVX  STALSX  RAMLSX 


There  are  8  linearly  independent  equations,  which  are  indicated  by  a 
line  joining  the  two  benchmarks,  one  in  each  region,  i.e.  P and  Qj 
between  which  (2.1)  is  formed.  Here,  WETZEL  and  PLTTVX  have  been  selected 
as  SMA  and  BMB,  again  rather  arbitrarily  but.  for  the  consideration  of 
bring  centrally  placed  in  each  region  (Figures  6.1  and  6.4)  to  give  op¬ 
timum  lengtns  for  leveling  from  BMA  or  BMB  to  cat;  centers .  The  .enter  i 
.cation  of  BMA  and  BMB  is  ..envenient,  as  we  nave  to  cons  i  be*  i  a  • 

.differences  W  (BMA.P.)  and  7W.  ( BMB , Q 1 )  by  leveling  to  enable 
f  the  8  equations  (2.1)  being  written  in  terms  cf  N(8MA.3MB.. 


o-  d  t. 


From  Figure  6.1,  it  is  reasonable  to  assume  that  AW, (BMA.P. )  may 
be  obtained  by  independent  routes  from  BMA  to  each  cap  center  P.  ,  and 
thus  V(-AW,)^  would  be  a  diagonal  matrix.  From  Figure  6.4,  if  we  as¬ 
sume  that  the  route  for  leveling  from  PLTTVX  to  RAMLSX  is  through  STALSX, 
then  V(eA W;)g  would  be  as  in  (2.13). 

The  oarticulars  for  the  Europe  -  USA  vertical  datum  connection  for 
-  =  3°,  2°,  1°  are  given  in  Table  6.1.  The  accuracy  estimate  •  ",W(BMA  ,BMB) 
is  also  given  for  the  three  cases.  To  highlight  the  effect  of  differing 
number  of  equations  for  the  three  values  of  ,  ,  rAW  is  quoted  without 
considering  any  geocentric ity  errors  in  SL4  coordinate  system,  and  with¬ 
out  considering  any  station  position  errors  or  any  leveling  errors.  The 
effect  of  these  later  errors  is  examined  in  Section  6.2.  The  normal 
gravitational  field  is  described  by  December  1981  potential  coefficients 
(Rapp,  1981).  The  standard  deviations  of  the  potential  coefficients  is 
considered  in  determining  rc'T  (see  equation  (3.2)).  The  standard  de¬ 
viation,  ceAg*  ,  of  modified  gravity  anomalies,  is  considered  as  2  mgals. 
The  modified  gravity  anomalies  were  considered  in  a  uniform  pattern  (see 
below  (4.7)*)  at  a  density  Ac  *  10'. 

We  note  from  Table  6.1  that  because  of  qreater  number  of  caps  being 
available  for  establishing  the  vertical  datum  connection,  anomaly  data 
in  capsize  of  .  =  1°  gives  :AW(BMA,BMB)  of  0.45  kgal .m,  which  is 
slightly  worse  than  0.37  kgal .m  for  ,.  =  2  .  But,  considering  that 
AW,(P.,P')  observations  by  leveling  are  required  from  cap  center  to  much 
lesser  number  of  anomalies  in  each  cao  (of  the  order  of  150  anom-il  ies  fv' 

-■  =  1°  to  the  order  of  500  for  ,  =2'),  capsize  of  .  =  1  appears  be 
adequate  for  establishing  the  vertical  datum  connection.  In  fact,  ’f 
we  use  =  0.5°  and  even  if  we  do  not  have  additional  caps,  we  get 
.'AW(BMA  ,BMB)  =  0.54  kgal  .m,  which  is  much  better  than  the  present,  likely 
accuracy  of  1.41  kgal .m  for  connection  between  coasta1  benchmarks  in 
the  two  regions  if  we  assume  the  r.m.s.  value  of  Sea  surface  tonocirap°v 
to  be  about  1  meter  . 

We  now  examine  the  effect  on  the  Europe  -  USA  vertica'  da*um 
nection  of  other  errors,  which  were  not  considered  in  Table  r . I . 
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Europe-  USA  Vertical  Datum  Connection  for  Various  Capsizes  . 
Normal  Gravitational  Field,  December  1981  Potential  Coefficients. 
Modified  Anomalies  Accuracy  o> An*  =  2mgals.  Anomaly  Density  10'. 
No  other  errors  considered. 
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for  P,:  (1)  WETZEL; (2)  GRASSE;  (3)  K00LAS;  (4)  F1NLAS 

for  Qj;  (1)  PLTTVX;  (2)  VERNLX;  (3)  BERLKX;  (4)  QU1NCX;  (5)  0VR79X;  (6)  GLD79X;  (7'  SANDIX 

(8)  FTDAVX;  ( 9 )  HOPLAS;  (10)  FLAGSX;  (11)  TLRBCX;  (12)  STALSX;  (13)  HAYSTX;  (14)  RAMLSX 

.1  for  detailed  information  on  these  SLR  stations. 


6.2  Variations  Due  to  Coordinate  System,  Position  and  Leveljna  Errors 

If  the  origin  of  the  SL4  coordinate  system  does  not  coincide 
with  the  center  of  mass  (geocenter)  of  the  earth,  it  results  in  errors 
in  the  positions  of  cap  centers,  whose  radial  component  is  given 
by  (2.9).  This  causes  additional  covariances  given  by  (2.10)  in 
the  variance  covariance  matrix  of  normal  grav i tational  potential 
differences  V(-AU)  of  (2.5).  The  effect  on  oaW(BMA ,BMB)  was  ex¬ 
amined  for  different  uncertainties  in  the  shifts  of  the  origin  of  SL4 
coordinate  system  from  the  geocenter,  as  character ized  by  the  standard 
deviation,  os  ,  of  the  shift.  ~s  was  varied  from  0  to  50  cm. 

Besides  the  shift  of  the  origin,  the  determination  of  station 
positions  have  errors  with  respect  to  the  SL4  coordinate  system.  The 
error  in  the  normal  gravitational  potential  is  primarily  caused  by  the 
radial  component  of  the  error  in  station  position,  and  is  given  by 
(2.7).  Assuming  no  correlation  between  radial  errors  of  different  cap 
centers,  the  variances  in  V(eAU)  are  given  by  (2.8).  The  effect  on 
cAW(BMA ,BMB)  was  examined  for  different  estimates  of  radial  errors,  as 
characterized  by  their  standard  deviation,  ocr.  :cr  was  also  varied 
from  0  to  50  cm. 

We  also  consider  the  effect  on  oAW(BMA,BMB)  of  leveling  errors 
in  estimating  potential  differences  AW,(BMA,P.)  and  AW?(BMB,Q.)  in 
relating  potential  differences  between  cap  pairs  P . ,  Q.  to  AW(BMA, 
BMB)  through  (2.1).  The  formation  of  V(cAAWj)  in  (2.5)  was  discussed 
in  Section  2.3  by  summation  of  corresponding  elements  in  V ( -  AW ; ) ^ 
and  V(-AW;)g  in  each  region.  The  latter  matrices  are  obtained  in 
the  form  of  (2.13)  depending  on  the  routes  of  leveling  connections 
from  BMA  (and  BMB)  to  cap  centers  P.  (and  Q ^ ) .  The  leveling  con¬ 
nections  for  the  Europe  -  USA  vertical  datum  connection  were  discussed 
in  Section  6.1,  and  specifically  shown  in  Table  6.1.  The  estimate  uf 
errors  in  AW  due  to  random  errors  in  leveling  is  aiven  by  (2.12): 

'  -  AW  =  k(0.l/T)  kgal.m,?,  in  10’  km 


The  effect  on  uAW(BMA  ,BMB)  was  examined  for  different  accuracy  es¬ 
timates  for  leveling  by  varying  k  from  0  to  4. 

The  accuracy  estimate  cAW(BMA  ,BMB )  for  the  Europe  -  USA  vertical 
datum  connection  is  tabulated  in  Table  6.2  for  the  variations  discussed 
in  this  section.  The  variation  in  oAW  is  tabulated  for  each  of  the 
three  errors  separately  for  the  cases  of  r  =  0.5“,  1°,  and  2°. 

We  note  from  Table  6.2  that  the  main  variation  in  cAW(BMA,BMB)  is 
due  to  the  anomaly  data  cap  size  as  discussed  in  Section  6.1.  The  ef¬ 
fect  of  uncertainties  in  the  shift  of  origin  (os)  of  the  SL4  system, 
ana  trie  effect  of  uncertainties  in  radial  position  (>:r),  on  :AW  is 
.-\..i'parat  iveiy  smaller.  The  effect  of  s  is  slightly  larger  as  com¬ 
pared  u 


effect  of  leveling  errors  on  :AW(BMA,BMB)  has  been  plotted  in 
t.5,  'nclcdir.g  some  additional  results  for  =  3°.  We  recall 

.  .i  *  ha  t  a  conservative  estimate  of  potential  differences 
■  .  :r'  ..ver  intra-continental  distances  is  given  by  0.1  ^ 

•  .o  .  ■  -f  '.T  xm  (see  (2.12)).  Even  if  this  estimate  is  doubled 

»  '  •  :a  ’  •’ .  *■  ne  effect  on  •AW  is  not  large  (see  last  3  rows  of 
'd-'-  n.  s  neans  that  even  very  conservative  estimates  of 

’  e  1  w  •;  errors  do  nut  significantly  affect  the  accuracy  estimates  of 
c^rc.ue  -  JSA  y-rt  id*.  datum  connection. 


Table  6.2  Variation  in  Accuray  Estimate  '.'V(BMA  ,8MB)  for  Europe-USA 
Vertical  Datum  Connection  due  to:  (a)  s.d.  (-s)  of  shift 
of  origin  of  SLR  stations  from  Geocenter;  (b)  s.d.  (  ■  r) 
of  radial  position  error  in  SLR  stations;  and  (c)  accuracy 
estimate  (ceAW0)  of  potential  difference  by  level ina  between 


cap  centers: 


Other  particulars  as  in  Table  6.1. 


=  k(0.1  vT)  kgal  .m, l  in  103  km. 


Capsize  v 


*  Caps  in  Europe 

4 

|  4 

4 

#  Caps  in  USA 

10 

- 

14 

| 

WETZEL 

WETZEL 

WETZEL 

1  PLTTVX 

|  PLTTVX 

PLTTVX 

Shift  i  Position 
as  (cm)  !  >:r  (cm) 


Level ing 
k 


Accuracy  Estimate  oaW(BMA,BMB)  kgal 


6.3  Optimum  Choice  of  Anomaly  Cap  Radius 


We  have  already  examined  the  effect  of  variation  in  anomaly  data 
cap  size  &  on  c<iW(BMA,BMB)  in  Section  6.1,  and  it  appears  that 
i  =  1°  or  P  =  0.5°  may  be  satisfactory  considering  the  comparatively 
slow  variation  in  oAW  due  to  other  errors  examined  in  Section  6.2. 

To  determine  an  optimum  ^ ,  we  need  to  adopt  conservative  estimates  of 
uncertainties  of  errors  in  Section  6.2;  and  consider  if  a  lower  value 
of  o  ,  requiring  lesser  field  work  of  establishing  AWJ,(P-j,P')  for 
modified  gravity  anomalies  (see  (4.9)*)  in  each  cap,  would  not  cause 
too  large  an  increase  in  oAW(BMA,BMB) . 

A  comparison  of  the  origin  of  SLR  station  coordinates  in  1980 
with  the  origin  determined  in  a  previous  solution  in  1979  gave  a  shift 
of  about  15  cm  (Smith  et  al . ,  1982).  Hothem  et  al .  (1982)  estimated 
the  relationship  of  co-located  SLR  station  coordinates  and  the  Doppler 
coordinate  systems.  After  allowing  for  about  4  meter  shift  in  the  Z 
direction  for  the  Doppler  coordinate  system,  they  estimated  the  dif¬ 
ference  in  the  origins  of  the  two  coordinate  systems  to  be  about  40 
cm  in  X,  70  cm  in  Y  and  10  cm  in  Z  directions.  Though  the  above  com¬ 
parison  was  strictly  not  with  the  SL4  solution  (Smith  et  al . ,  1982),  and 
the  origin  of  the  SL4  system  may  be  located  more  precisely  with  the 
geocenter,  we  may  assume  a  conservative  estimate  of  os  in  (2.11)  as 
30  cm. 

The  formal  standard  deviation  of  the  annual  mean  SLR  station 
heights  in  the  SL4  system,  based  on  the  random  error  estimate  of  laser 
ranges,  is  of  the  orderof  1  to  2cm(D.  Chri stodoul idi s ,  personal  com¬ 
munication,  July  1982).  The  standard  deviation  of  the  monthly  solu¬ 
tions  of  SLR  station  heights  from  Lageos  data  in  1979-80  varies  from 
7  to  12  cm  (Smith  et  al . ,  1982).  A  conservative  estimate  of  radial 
errors  in  SL4  station  position  may  be  assumed  as  oer  =  10  cm. 

A  conservative  estimate  cf  errors  in  potential  difference  through 
leveling  over  intra-continental  distance  was  discussed  in  (2.12)  as 
oeAW  =  0.1/T  kgal.m,  i  in  101  km.  We  give  in  Table  6.3  the  accuracy 
estimate  of  Europe  -  USA  vertical  datum  connection  for  '•  'W  as 


0.1  ,7,  also  for  0.2  71;  and  to  emphasize  that  even  very  large  leveling 
errors  do  rot  significantly  affect  ;AW(BMA ,BMB) ,  we  also  give  results 
for  .:eAW£  as  0.4  kqal.m,  t  .in  103  km.  Anomaly  data  in  a  uniform 
pattern  around  each  cap  center  is  considered  for  y  =  1  ’  and  0.5';  and 
also  given  for  comparison  for  o  =  2C  as  we1!  as  y  =  0.25'  and  .  =  0.17 
The  anomaly  density  Ay  was  kept  at  10',  as  a  higher  density  does  not 
cause  any  significant  improvement  in  cel  (Figure  4.2).  However,  to 
ensure  that  at  least  3  ring  averages  are  used  for  estimating  f(P.) 
or  T(Q  )  ,  the  density  was  increased  for  v  =  0.25°  and  y  =  0.1°.  The 
number  of  modified  anomalies  iri  a  uniform  pattern,  and  :cT  at  each  cap 
center  are  also  given  in  Table  6.3. 


The  accuracy  estimates  ;7W(BMA ,BM3 )  in  the  last  three  rows  of 
Table  6.3,  and  additional  values  for  ,  =  3',  are  shown  in  Figure  6.6. 

We  note  that  a  conservative  accuracy  estimate  for  Europe  -  USA  vertical 
connection  with  anomaly  data  caps  of  spherical  radius  of  V  and  0.5' 
is  0.50  and  0.58  kqal.m,  which  changes  to  .52  and  .59  kgal .m  respectively- 
even  if  the  estimate  of  leveling  errors  is  doubled.  The  accuracy  es¬ 
timates  will  be  further  improved,  when  additional  laser  ranging  stations 
in  Europe  in  Spain,  Switzerland  and  Greece  are  included  in  the  lageos 
SLR  network. 


Considering  that  the  requirement  of  field  work  in  each  cap  for 
providing  leveling  connections  from  cap  center  to  the  modified  anomalies 
(see  (4.9)*)  is  very  considerably  reduced  for  ,  =  0.5'  as  compared  to 
=  1'  (by  about  one-third  the  number  of  anomalies  over  one- fourth  the 
area),  the  optimum  choice  for  anomaly  data  cap  radius  appears  to  be 
.  -  0.5'.  This  leads  to  a  very  modest  data  requirement  of  about  50 
to  60  anomalies  at  a  rough  spacing  of  15  to  20  km  around  the  SLR  sfatirnc 
in  a  diameter  of  about  120  km.  Such  anomalies,  with  level inq  connection 
to  tne  cap  center,  are  already  available  around  most  SLR  stations,  or 
may  oe  established  with  modest  field  work.  These  anomalies  ■>  i  .  then 
be  used  to  predict  the  modified  anomalies  in  a  uniform  pattern,  oe- 
Sv.r'bed  below  (4.7)*. 
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Table  6.3  Accuracy  Estimate  cAW(BMA,BMB)  for  Europe-USA  Vertical 
Datum  Connection. 

Optimum  Choice  of  Anomaly  Data  Capsize 
Normal  Gravitational  Field,  December  1981  Potential  Coefficients, 
s.d.,  os,  of  shift  of  origin  of  SLR  stations  from  geocenter 
o.  =  30  cm 

s.d,  oer,  of  radial  position  error  in  SLR  stations:  •  r  -  10  cm 
s.d.,  oeAW-,  of  potential  difference  errors  by  leveling  between 

cap  centers:  ocAW  =  k (0 . 1  vT)  kgal  .m,  ■’  in  103  km,  k=l,  2,  4. 
s.d,  oeAg*  =  2  mgals. 


Capsize  u>° 

P.:  4  Caps  in  Europe 
Q-:  #  Caps  in  USA 


Anomaly  density ,Ai|j ,  in 
uniform  pattern 

Approx,  diameter  of  anomaly 
data  cap  (km) 

fr  Anomalies  per  cap 

#  Anomalies  ring  averages 
per  cap 

Accuracy  of  predicted  T  in 
a  cap  (j?T  -  kgal.m) 


1°  i  0. 
4  j  4 


WETZEL 

PLTTVX 

WETZEL 

PLTTVX 

WETZE 

PLTTi 

10' 

7.5' 

3 

110 

56 

22 

43 

19 

19 

4 

3 

3 

.82 

.95 

1 

Potential  Difference  Errors 
Through  Leveling 
>cAW5(BMA,P.);  ocAW„  (BMB ,Q  • ) 


O.lv^Tkgal.m,  l  in  103km 


oAW(BMA ,BMB)  kgal.m 


ACCURACY  aAW  (BMA.BM  B)  kgal. 


Figure  6.6:  Accuracy  Estimate  oW(BMA,BH8)  of  Europe  -  USA  Vertical 

Datum  Connection. Optimum  Choice  of  Anomaly  Cap  Radius  ;°, 
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ACCURACY  ESTIMATE  FOR  OTHER  VERTICAL  DATUM  CONNECTIONS 


We  apply  in  this  section  the  conclusions  of  Section  6.3  to  determine 
the  accuracy  estimate  for  other  vertical  datum  connections.  The  datum 
connections  considered  were  that  of  U.S.  respectively  with  Australia, 

Hawaii,  Peru  (in  Sou*h  America),  Bahamas  and  Bermuda  (in  mid-Atlantic 
off  the  U.S.  East  Coast).  There  are  two  SLR  stations  at  Orroral  and 
Yarragadee  on  the  East  and  West  coast  of  Australia,  see  Figure  3.1;  and 
one  SLR  station  each  at  other  sites,  assuming  that  a  satisfactory  intra¬ 
continental  leveling  connection  may  not  be  available  between  the  two 
SLR  station  sites  in  South  America  at  Arequipa,  Peru  and  Natal,  Brazil. 

In  all  cases,  14  SLR  station  sites  were  considered  in  U.S.  First,  the 
accuracy  estimates  oAW(BMA ,BMB)  were  obtained  with  anomaly  data  cap 
radius  V  around  each  SLR  station  as  0.5°;  and  as  'AW(BMA  ,BMB)  was 
comparatively  larger  than  in  Section  6  due  to  smaller  number  of  caps  in 
the  vertical  datum  connection,  oAW(BMA,BMB)  was  also  obtained  with 
V  =  lc.  In  both  cases  of  '*>  =  0.5°  and  P  =  1°,  cAW(BMA  ,BMB)  was  obtained 
with  different  variations  cs,  cer,  c-:AW  as  in  Section  6.2  (Table  6.2). 

We  quote  below  in  Table  7.1  only  the  results  for  the  standard  deviations 
as  ,  shift  of  SLR  coordinate  system  from  geocenter,  as  30  cm;  -rr,  radial 
position  error  of  SLR  stations,  as  10  cm;  and  -cAW,,  ,  potential  dif¬ 
ference  between  cap  centers,  as  0.1  VT  legal  .m,  where  •  is  the  length 
of  level  line  in  103km.  jAW(BMA,BMB)  is  also  given,  when  none  of  the 
above  errors  are  considered,  i.e.,  each  of  os,  ~r,  :.AW„  set  as  zero. 

This  allows  for  a  ready  comparison  with  results  in  Section  6,  where  we 
had  4  SLR  sites  in  Europe,  while  we  have  now  only  2  SLR  sites  for  the 
vertical  datum  connection  with  Australia,  and  only  1  SLR  site  for  other 
vertical  datum  connections.  The  accuracy  estimate  AW(BMA,BMB)  were 
almost  similar  for  the  latter  4  datum  connections,  i.e.  of  U.S.  respectively 
with  Hawaii,  Peru,  Bahamas,  Bermuda.  Hence  the  results  for  only  one  of 
these,  i.e.  U.S.  -  Bermuda  vertical  datum  connection  have  been  given  in 
Table  7.1. 


laDle  7.1  Accuracy  Estimate  of  j .'1 . and  Other  Vertical  Datum  Connections. 

•s  =  s.d.  of  shift  of  origin  of  SLR  stations  from  geocenter;  =  30  civ 
•cr  =  s.d.  of  radial  position  errors  in  SLR  stations;  v-;r  =  10  cm 
cczlW  .  =  s.d.  of  potential  difference  errors  by  leveling  between  cap 
centers  in  a  region 

=  O.lvT  kgal  .m,  t  in  10 3  km. 

Other  particulars  as  in  Table  6.3 


Vertical  Datum  Connections 


No.  or  caps  in  U.S. 

No.  of  Caps  in  Other  Datum 


I  U.S. -Europe 


U.S. -Austral ia 


U.S. -Bermuda 


Anomaly  Data  Capsize 


Accuracy  Estimate  cAW(BMA  ,BMB)  kgal  .r 


and  cclW,  =  O.lvT  kgal.  m 
l  in  TO3  km 


Accuracy  Estimate  oaW(BMA,BMB)  kgal  .r 


Anomaly  Data  Capsize  o  =  r 

vs  *  r  =  -t-AWj,  =  0 

s  =  30  cm,  r  =  10  cm, 

and  =  0.1  v'Tkgal.m, 

tin  lO3  km 


*  Similar  results  were  obtained  for  other  vertical  datum  connections,  with 
only  one  SLR  station  in  the  second  region,  i.e.  U.S.  and  respectively 
Bahamas,  Hawaii,  Peru. 


We  note  from  Table  7.1  that  the  accuracy  estimate  of  the  vertical 
datum  connection  is  poor  (.73  kqal.ni)  even  with  .  =  l'\  when  we  have 
only  one  SLR  station  in  the  second  region.  When  the  number  of  SLR 
stations  in  the  second  region  increases  to  2  and  4  respectively,  the 
accuracy  of  the  vertical  datum  connection  is  0.62  kgal.m  with  .  =  1  , 
and  0.58  kgal.m  with  0  =  0.5°  respectively.  The  number  of  stations 
available  for  the  vertical  datum  connection  has  a  greater  effect  for 
improving  the  accuracy  of  datum  connection's  compared  to  the  anomaly 
data  cap  size  . 

To  confirm  the  above  conclusion,  we  show  in  Table  7.2  the  accuracy 
estimate  of  U.S. -Europe  vertical  datum  connection  with  =  0.5°  for  three 
cases,  when  the  number  of  stations  in  Europe  are  taken  as  2,4  and  7 
respectively.  For  the  first  case,  we  consider  only  two  SLR  stations  in 
Europe  at  Wetzel! ,  Germany  and  Grasse,  France  (see  note  below  Table  6.1 
and  Figure  6.1).  For  the  second  case,  we  consider  4  stations  in  Europe 
as  in  Section  6.  For  the  third  case,  we  consider  three  additional 
stations  in  Europe  at  Dionysos,  Greece;  Zimmerwald,  Switzerland;  and 
San  Fernando,  Spain.  The  coordinates  of  these  additional  stations  are 
not  currently  available  in  the  SL4  (Smith  et  al.,  1982)  solution  of 
Lageos  SLR  network,  but  these  coordinates  are  likely  to  be  available  in 
subsequent  solutions. 

Table  7.2  Accuracy  Estimate  of  U.S. -Europe 

Vertical  Datum  Connection.  Variation 
Due  to  Number  of  SLR  stations  in  Europe. 

Anomaly  Data  Capsize  y  =  0.5° 

os  =  30cm,  otr  =  10  cm,  •:  AW,  =0.1  kgal.m 

(see  explanation  in  Table  7.1) 

Other  particulars  as  in  Table  6.3. 


No.  of  Caps  in  U.S. 
No.  of  Caps  in  Europe 


14 

2 


Anomaly  Data  Capsize  ;  =  0.5  ; 

s  =  30  cm,  r  =  10  cm,  and 
'-■'-W-  =0.1  /f  kgal.m,  ■  in  10  ’  kn ' 


Accuracy  Estimate 


AW(BMA.BMB)  kga 
.  74 


.49 


.58 


A  comparison  of  the  last  row  and  the  first  columns  of  Tables  7.1  and 
7.2  shows  that  the  accuracy  estimate  of  U.S. -Europe  vertical  datum  is  about 
the  same  when  we  use  .  =  0.5 5  but  have  three  additional  stations  in  Europe 
as  compared  to  =  1  but  without  the  three  additional  stations.  However, 

there  is  a  large  increase  in  field  work  of  leveling  from  the  cap  centers 
to  the  anomalies  in  the  latter  case  of  v  -  1'-  Hence,  the  number  of 
stations  available  for  the  vertical  datum  connection  has  a  greater  effect 
for  improving  the  accuracy  of  datum  connection,  as  compared  to  the  anomaly 
data  capsize 
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SUMMARY  AND  CONCLUSIONS 


We  have  investigated  the  procedures  recommended  by  Colombo  (1980a)  for 
establishing  a  'world  vertical  network',  and  have  applied  the  results  to 
determine  the  current  accuracy  estimate,  in  particular,  of  the  Western 
Europe- USA  vertical  datum  connection. 

The  principle  of  the  method  is  to  determine  the  potential  difference 
between  two  'bench-marks'  BM  A  and  BM  B  in  the  two  regions  A  and  B,  without 
a  precise  knowledge  of  the  mean  sea  level,  i.e.  without  a  precise  knowledge 
of  the  'sea  surface  topography'  or  the  geoid,  in  particular  the  zero  degree 
term  in  the  earth's  gravitational  potential,  the  effect  of  which  nearly 
cancels  out  in  determining  the  potential  difference.  The  estimate  of  the 
potential  difference  AW(BMA,BMB)  is  more  precisely  obtained  through  usual 
least  squares  technique  by  redundant  observations,  when  we  consider  the  po¬ 
tential  difference  between  several  pairs  of  benchmarks  P .  and  Q ^ ,  respectively 
in  the  two  regions  A  and  B,  and  relate  the  potential  difference  of  each 
pair  to  AW(BMA,BMB)  by  determining  AW(BMA,  P..)  and  £W(BMB,  Qj)  through 
classical  leveling  and  gravity  observations. 

Hence,  the  requirements  for  establishing  a  vertical  datum  connection  are: 

(a)  Precise  station  positions,  particularly  in  the  radial  component,  in  an 
earth-centered  coordinate  system,  to  determine  the  normal  gravitational  poten¬ 
tial  U(P  - )  or  U(Q.)  and  the  centrifugal  potential  T(p  )  or  !((}•)--  see 

1  J  *  J 

(2.1).  The  possible  shift  of  the  origin  of  coordinate  system  from  the  geo¬ 
center  is  considered  through  correlated  errors  in  the  normal  gravitational 
potential  estimates  of  the  stations  (see  (2.10)).  The  random  errors  in 
station  positions  in  the  established  coordinate  system  are  considered  through 
(2.7)  and  (2.8).  Any  systematic  effect  of  a  scale  factor,  the  radial  com¬ 
ponent  of  all  station  positions  being  larger  or  smaller,  is  nearly  cancelled 
out  in  determinina  Dotential  differences. 


(b)  An  accurate  description  of  the  earth's  gravitational  field  in  terms  of 
potential  coefficients  in  a  spherical  harmonic  expansion  to  a  high  degree 
N.  Th^  enab1  es  comouta tic"  of  U  through  (2.6).  The  errors  in  the 
Potential  coefficients  a-e  considered  in  defining  the  anomalous  potential 


T  in  (3.2)  through  the  estimated  standard  deviations  cf  the  potential 
coefficients  in  (3.4)  to  (3.7).  The  potential  coefficients  f,r  degrees 
higher  man  N  are  modelled  through  anomaly  degree  variance  models--  see 

(3.8)  to  (3.11). 

(c)  Gravity  anomalies,  residual  to  degree  N  ,  at  an  optimum  density  and 
extent,  or  capsize,  around  each  station  to  predict  T'P.)  or  T(Q^)  through 
(4.1).  The  available  anomalies  are  first  used  to  predict  a  set  of  anomalies 
in  a  uniform  pattern  (see  below  (4.7)*)  for  an  efficient  algorithm  util¬ 
izing  'ring  averages'  (4.4),  which  greatly  reduces  the  size  of  matrices 

in  the  collocation  formulas  (4.1).  The  choice  of  optimum  density  and  ex¬ 
tent  of  the  anomalies  in  the  pattern  are  important  for  a  precise  estimate  of 
T  (Section  4.3).  Tne  propagation  of  random  errors  in  the  anomalies  to 
accuracy  estimate  of  T  is  considered  in  Section  4.5.  A  systematic  error 
in  gravity  anomalies  is  nearly  cancelled  in  tne  potential  differences  as 
the  same  density,  extent,  and  pattern  of  anomalies  is  used  around  each 
station . 

(d)  Estimate  of  potential  differences  through  leveling  between  the  locations 
of  gravity  anomalies  in  the  'cap'  and  the  station  or  cap  center.  This  en¬ 
ables  the  computation  of  'modified 'grav ity  anomalies  (Section  4.2),  such 
that  the  modified  anomalies  are  referenced  to  the  normal  potential  of  the 
cap  center  instead  of  being  referenced  to  the  potential  of  the  geoid,  see 

(4.8) *,  (4.9)*.  The  precise  knowledge  of  the  geoid  is  therefore  not  re¬ 
quired.  In  subpara,  (c)  above  ,  modi fied  gravity  anomalies  are  used  for 
predicting  T  ,  but  as  the  covariances  in  tne  collocation  formula  are  com¬ 
puted  utilizing  anomaly  degree  variance  model,  the  accuracy  estimate  of 

T  are  corrected  as  in  (4.13). 

(e;  Estimate  of  potential  differences  pW  ,  (3MA  .P.. )  and  SW  ,  (8MB ,Q  ■ ) ,  through 
lev->I  ing  and  rravity  observations,  over  intra- continental  distances  in  ea.h 
r-‘<!  on  A  d’’H  8.  The  effect  of  random  errors  in  a*-e  c or, i •  r  !  •» n 

‘  ■: ,  ‘  •  ri  2.3.  Any  systematic  error  in  '.W  .  which  is  dependent  :i  >  ;  >vu 
u • e  0  t  re  is  (BMA.P.)  or  (8MB, Q.),  will  become  apparent  when  equations  (2.11 
j;  -i)  are  f r.rmed  . 


(f)  An  adequate  number  of  stations  for  establishing  the  vertical  connection 
to  enable  a  more  precise  estimate  of  AW(BMA,BMB).  The  number  of  equations 
(2.1)  or  (2.3)  are  one  less  than  the  total  number  of  stations  in  both  re¬ 
gions  A  and  B. 

8.1  Recommendations  of  Present  Study 

The  above  requirements  are  already  fulfilled  to  a  sufficient  extent 
to  attempt  the  vertical  datum  connection  between  Western  Europe  and  USA. 

A  modest  amount  of  field  work  may  be  necessary  before  actual  computations 
are  done.  We  first  list  the  current  possibilities,  and  identify  the  er¬ 
rors  whose  effects  on  the  vertical  datum  connection  was  investigated  in 
this  study;  and  also  briefly  mention  the  highly  near-future  improvements: 

(a)  If  we  count  the  repeated  occupations  of  adjacent  marks  in  a  vicinity 
as  one  station,  i.e.  using  only  one  of  the  marks,  we  have  29  station  posi¬ 
tions  (Figure  3.1)  in  the  Lageos  Satellite  ranging  network  in  the  SLA  co¬ 
ordinate  system  (Smith  et  al.,  1982).  We  used  a  very  conservative  estimate 
of  the  standard  deviation  of  the  shift  of  SL4  coordinate  system  from  the 
geocenter  as  30  cm,  and  we  used  the  standard  deviation  of  the  radial  com¬ 
ponent  in  the  station  positions  in  the  SL4  coordinate  system  as  10  cm 

(see  Section  6.3).  These  are  likely  to  be  substantially  improved  in  the 
subsequent  SL5  solution.  However,  this  is  not  a  critical  factor  and  may 
improve  the  accuracy  estimate  of  the  vertical  connection  by  only  a  marginal 
amount,  say  .03  kgal .m  (see  Table  6.2). 

(b)  We  used  the  earth's  gravitational  field  model  developed  by  Rapp  (1981), 
which  was  termed  in  this  study  as  December  1981  potential  coefficients  field, 
and  was  complete  up  to  degree  180.  The  standard  deviations  of  the  potential 
coefficients  are  perhaps  slightly  pessimistic  (Lachapelle  and  Rapp,  1982, 

p.  4).  Any  improvement  in  the  accuracy  of  the  grav itational  field  will  re¬ 
sult  in  significant  improvement  of  the  accuracy  estimate  of  the  vertical  con¬ 
nection--  see  Figure  4.3  and  Table  4.3.  A  revised  estimate  of  the  gravity 
field,  with  accurately  determined  low  degree  (say  N  10)  potential  coef- 
ficients^may  improve  the  accuracy  estimate  of  the  vertical  connection  by 
about  0. 1  kgal .m. 
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(c)  The  optimum  density  of  gravity  anomalies  is  at  spacing  of  about  10', 
roughly  15  to  20  km  apart -- see  Figure  4.2  and  Table  4.2.  The  optimum  ex¬ 
tent  of  anomalies,  or  the  radius  c  of  the  anomaly  data  cap  around  each 
station,  is  obtained  by  a  judicious  balance  of  the  practical  feasibility 
of  the  increased  field  effort  with  larger  -v'  versus  the  lowered  precision 
of  the  accuracy  estimate  of  the  vertical  connection  with  smaller  ..  An 
important  conclusion  is  that,  unlike  Colombo's  original  proposal  of  .  =  5:, 

*  may  certainly  be  reduced  to  1°  (Sections  6.1  and  6.2)  without  significant 
dete r i ora t i on  in  vertical  connection  accuracy,  particularly  because  more 
stations  are  available  with  reduced  cap  size  (compare  Figures  6.2,  6.3  and 
6. A).  A  further  reduction  in  y  to  0.5°  is  discussed  below.  A  very  en¬ 
couraging  conclusion  is  the  lack  of  sensitivity  of  the  vertical  connection 
accuracy  to  the  random  errors  in  gravity  anomalies  (Figure  4.5  and  Table 
4.4).  With  several  SLR  stations  located  in  the  coastal  areas,  and  anomaly 
data  cap  covering  marine  areas,  it  is  very  helpful  to  know  that  even  large 
random  errors  of  10  mgals  (1-j)  in  anomalies  will  lower  the  accuracy  of  ver¬ 
tical  connection  only  by  about  0.01  kgal.m. 

(d)  Gravity  observations  at  approximate  spacing  of  15  to  20  km  may  be 
available  around  most  stations,  but  leveling  observations  may  be  required 
to  determine  potential  differences  from  the  station  at  cap  center  to  the 
locations  of  anomalies.  When  we  consider  that  in  case  of  ^  =  0.5',  about, 
one- third  anomalies  are  required  to  be  connected  by  leveling  over  one- fourth 
the  area  around  each  station  (Table  6.3),  as  compared  to  the  case  of  , .  =  l  , 
while  the  accuracy  of  vertical  connection  is  lowered  by  only  about  0.07  kaal.m. 
it  appears  that  the  optimum  extent  of  gravity  anomalies  is  a  circle  of  di¬ 
ameter  of  about  120  km  corresponding  to  4,  =  0.5’. 

(e)  The  vertical  connection  accuracy  is  quite  insensitive  to  even  very 
conservative  estimates  of  random  errors  in  leveling  over  intra-continental 
distances,  see  Figures  6.5  and  6.6  and  equation  (2.12).  The  currently  avail¬ 
able  estimates  of  AW,,  (BMA ,  P  )  and  AW,(BMB,Q,)  are  therefore  adequate 

*-  1  ■  j 

ter  the  vertical  datum  connections.  This  is  a  very  encouraging  conclusion 
in  view  of  concern  over  accuracy  of  leveling  ever  intra-continenta i  distances, 
e.g.  Kumar  and  Soler  (1981).  Further,  the  effect  of  systematic  ei rors  in 

•  should  become  apparent,  by  the  examination  of  residuals  in  ennations  '?.  >  1 . 
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(f)  The  vertical  connection  accuracy  improves  sharply  with  the  increase  in 
the  total  number  of  stations,  see  the  last  two  rows  of  Table  6.1.  We  have 
presently  used  4  stations  in  Europe  and  14  stations  in  USA.  This  should  be 
increased  in  the  very  near  future  by  another  three  stations  when  the  co¬ 
ordinates  of  laser  ranging  stations  in  Spain,  Switzerland  and  Greece  are  also 
solved  for  in  the  Lageos  network.  The  number  of  stations  would  increase 
further  in  the  next  couple  of  years  with  the  utilization  of  Transportable 
laser  ranging  stations  (NASA,  1982). 

The  present  accuracy  estimate,  with  the  data  used  in  this  study  (see 

(a)  to  (f)  above),  for  the  vertical  datum  connection  between  Western  Europe 
and  USA  is  about  60  kgal.cm,  see  Table  6.3.  The  improvement  in  this  es¬ 
timate  would  come  primarily  from  increased  number  of  stations  and  more  ac¬ 
curate  estimates  of  earth's  gravity  field,  subparas  (f)  and  (b)  above. 

The  very  near  future  prospects  for  the  vertical  datum  connection  accuracy 
are  likely  to  be  about  50  kgal.cm.  This  is  about  one-third  of  the  classical 
vertical  connection  accuracy  through  tide  gage  determinations  due  to  sea 
surface  topography. 

The  various  requirements  detailed  in  Section  8  are  already  available, 
or  easily  accessible.  The  main  additional  field  effort  is  the  estimate  of 
potential  differences  at  each  SLR  station  by  establishing  leveling  con¬ 
nections  from  the  station  to  the  locations  of  about  50  gravity  anomalies  at 
a  spacing  of  15-20  km  in  a  circle  of  diameter  of  about  120  km  around  the 
station. 

8_. 2  Recommendations  for  Further  Studies 

The  following  topics  were  not  covered  in  the  present  investigation,  and 
should  be  studied  further: 

(a)  An  optimum  algorithm  for  predicting  anomalies  in  a  uniform  pattern, 
see  below  (4.7)*,  around  each  SLR  station  from  the  nearest  available  anom¬ 
alies. 

(b)  An  optimum  procedure  for  computing  modified  gravity  anomalies  from  the 
available  gravity  observations;  and  consideration  of  the  corrections  to  be 
applied  to  such  observations. 

(c)  The  problem  of  determining  the  absolute  value  of  geopotentia1  at  dM  A 
or  8^  8,  with  better  or  comparable  accuracy  to  the  potential  diffe»>?rce 

W(3VA,BK3)  of  the  vertical  datum  connection. 


(d)  The  possibility  and  the  accuracy  of  sea  surface  topography  estimates 
by  connecting  the  coastal  SLR  stations  to  tide  gages  and  comparing  with  the 
mean  sea  level  determinations  at  those  tide  gages. 

A  better  insight  into  the  establishment  of  the  'World  Vertical  Network' 
will  be  obtained  by  initiating  the  computations  of  the  vertical  datum 
connection  between  Western  Europe  and  USA.  Most  of  the  data  is  easily 
accessible,  and  a  very  modest  field  effort  is  needed  to  obtain  the  re¬ 
maining  data.  A  very  intriguing  question  is  the  possibility  of  revealing 
any  systematic  errors  in  leveling  observations  over  intra-continental  dis¬ 
tances  by  examining  the  residuals  of  equations  (2.1).  This  may  be  the 
appropriate  time  to  start  such  investigations  in  view  of  the  proposed  re¬ 
definition  of  the  vertical  networks  in  North  America. 
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